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As  described  in  previous  reports,  this  AASERT  funding  has  provided  support  for  graduate 
student  training  in  high-speed  optical  flow  diagnostics.  Funds  were  provided  for  the  purchase  of 
a  20  MHz  rate  framing  camera  (Hadland  Imacon  792),  which  was  used  to  capture  a  pair  of 
rapidly  sequenced  Rayleigh  scattering  images,  using  a  double-pulsed  Nd:YAG  laser  system 
purchased  form  parent  grant  AFOSR  support.  During  the  current  reporting  period,  work  has 
progressed  on  a  synergistic  program  to  develop  a  "burst  mode"  laser  system,  which  will, 
ultimately,  provide  the  capability  to  obtain  approximately  30  rapidly  sequenced  images.  This 
work  is  closely  coupled  with  the  parent  grant  (#F496200-92-J-0217)  and  a  collaborative  AFOSR 
program  with  Profs.  Alexander  Smits  and  Gany  Brown  (F49620-93- 1-0064). 

During  the  current  reporting  period,  AASERT  funds  have  been  used  to  support  a  new  Air  Force 
program  to  explore  the  potential  of  a  radiatively-driven  hypersonic  wind  tunnel  facility.  The 
AASERT  graduate  student  (Mr.  George  Williams)  has  been  responsible  for  one-dimensional 
modeling  of  non-isentropic  expansions  under  highly  non-ideal  gas  conditions.  His  work  is 
documented  in  detail  in  the  attached  two  reprints  (Appendix  A),  and  we  will,  therefore, 
summarize  the  principal  results  below. 

A.  Results:  Radiative  Driven  Wind  Tunnel 

A  convenient  way  of  viewing  the  operation  of  such  a  facility  is  with  the  help  of  a  Mollier 
diagram  (Fig.  1).  In  this  diagram,  the  enthalpy  is  plotted  versus  the  entropy.  For  a  conventional 
wind  tunnel  facility,  the  expansion  is  assumed  to  be  isentropic,  so  the  static  temperature  and 
pressure  and  the  Mach  number  in  the  test  chamber  mandate  a  unique  temperature  and  pressure 
for  the  gas  in  the  plenum.  For  example,  for  a  Mach  16  flow  with  a  static  temperature  of  200  K 
and  a  static  pressure  of  10-3  atm  (150,000’-point  C  on  the  diagram),  the  plenum  must  be  held  at 
a  temperature  of  6300  K,  and  at  a  pressure  of  40(X)  atm  (Point  A).  Clearly,  at  this  temperature, 
air  is  highly  dissociated  and  severe  problems  with  containment  and  throat  erosion  occur.  In 
contrast,  the  radiatively-driven  wind  tunnel  begins  with  a  plenum  pressure  on  the  order  of  30,000 
atm,  but  a  plenum  temperature  on  the  order  of  1000  K.  The  point  on  the  Mollier  diagram  (Point 
B)  associated  with  this  initial  condition  corresponds  to  significantly  lower  entropy  and  lower 
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enthalpy  than  that  of  the  isentropic  wind  tunnel.  That  remaining  enthalpy  and  entropy  is  then 
added  downstream  of  the  throat  using  the  radiative  source. 

George  Williams'  work  has  focused  on  development  and  application  of  a  one-dimensional,  real 
gas  equation-of-state,  flow  model  which  incorporates  energy  addition  and  uncoupled  NOx 
kinetics.  He  has  explored  a  variety  of  potential  energy  sources,  including  the  HF  laser  at  2.7 
microns,  which  couples  to  air  through  the  order  lOO's  ppm  of  naturally  occurring  CO2;  60  GHz 
microwave  radiation,  which  couples  directly  to  molecular  oxygen;  and  high  energy  (2-20  MeV) 
electron  beams,  which  couple  by  inelastic  electron-neutral  collisions.  In  the  case  of  laser  energy, 
the  beams  are  assumed  to  propagate  in  the  upstream  direction  from  an  initial  location  in  the 
vicinity  of  a  test  section.  The  model  incorporates  refraction  effects,  due  to  the  large  index-of- 
refraction  gradients  which  would  exist  in  a  thermal  boundary  layer,  and  which  cause  the  beams 
to  be  steered  toward  the  flow  centerline. 

A  detailed  description  of  the  individual  elements  of  the  ID  code,  as  well  as  a  series  of  illustrative 
results  for  a  full-scale  facility  using  HF  laser  radiation  has  been  recenUy  published  in  the  AIAA 
Journal  (reprint  attached).  In  addition  to  fluid  phenomena  (i.e.,  that  which  can  be  described 
through  continuity  equations),  the  radiatively-driven  wind  tunnel  concept  requires  detailed 
knowledge  of  a  variety  of  molecular  phenomena,  in  particular,  optical  absorption  and  energy 
thermalization.  While  the  phenomenology  of  these  effects  is  well  understood,  many  details  have 
not  been  explored  in  the  temperature  and  pressure  regimes  relevant  to  the  radiative  wind  tunnel 
concept.  We  have  begun  an  experimental  and  modeling  effort  to  address  the  most  important 
issues,  preliminary  results  of  which  were  presented  at  the  30th  AIAA  Thermophysics 
Conference,  June  1995  (San  Diego,  California).  (Reprint  attached.) 

RecenUy,  our  Air  Force  sponsors  have  expressed  interest  in  smaller  scale  demonstration 
experiments  which  would  serve  to  demonstrate  the  principle  and  to  validate  modeling  efforts 
(both  ID  and  axisymmetric).  These  experiments  would  require  sources  capable  of  output  power 
on  the  order  of  1  MWatt  for  time  durations  on  the  order  of  0.1- 1.0  seconds.  Possible  sources 
include  a  high  power  DF  laser  at  White  Sands,  New  Mexico,  a  2  MeV  electron  beam  facility  at 
Sandia  National  Laboratories  in  Albuquerque,  New  Mexico,  and  a  1  MWatt,  110  GHz 
microwave  gyration,  at  either  General  Atomics  (San  Diego,  CA),  or  Varian  (Palo  Alto, 
California). 
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Figures  2-5  shows  some  recent  modeling  predictions  for  a  2  MeV  electron  beam  source  in  which 
the  final  total  enthalpy  corresponds  to  approximately  200%  of  the  stagnation  enthalpy.  The 
stagnation  conditions  are  1000  atm  P,  500  K,  which  corresponds  to  an  enthalpy  of  approximately 
0.5  MJ/kgm  of  air. 

Figures  2-5  show  the  ID  code  results  for  nozzle  radius,  free  stream  velocity,  and  free  stream 
static  pressure  and  temperature,  respectively.  In  this  case,  the  code  was  constrained  to  solve  the 
continuity  equations  for  a  constant  Mach  number  expansion  until  a  static  temperature  of  400  K 
was  achieved.  At  this  point,  the  expansion  was  constrained  to  isothermal. 

It  can  be  seen  that  a  relatively  gentle  expansion,  from  diameter  2.4  to  10  mm  over  a  length  of 
0.10  m,  is  required  to  accommodate  the  increased  enthalpy  from  the  electron  beam. 

Figures  2-5  are  meant  to  serve  as  illustrative  examples  only.  In  particular,  we  are  now  in  the 

process  of  assembling  an  axisymmetric,  two-dimensional  code  which  will  be  the  next  generation 
design  tool. 


B.  Pulse  Burst  Laser 

While  there  has  been  enormous  progress  in  recent  years  in  the  development  and  application  of 
diagnostic  imaging  techniques,  the  ability  to  capture  time-evolving  or  volumetric  information  is 
severely  constrained  by  limitations  of  available  laser  technology.  As  a  general  rule,  in  the  gas 
phase,  signal  levels  for  demonstrated  imaging  diagnostic  techniques  are  sufficiently  small  that 
the  measurements  require  high-power,  nanosecond  duration  pulsed  sources,  such  as  Q-switched 
solid-state  (principally  Nd:  YAG)  or  excimer  lasers.  The  pulse  repetition  rate  of  these  devices  is 
limited  to  the  range  10-30  pulses/second  for  solid-sate  lasers,  and  10-300  pulses/second  for 
excimer  lasers.  The  current  state-of-the-art  is  represented  by  double-pulsed  Nd:YAG  lasers, 
which  may  be  used  to  produce  a  pair  of  excitation  sheets,  separated  in  time  between  1-80 
microseconds,  to  capture  two  sequential  images.  We  are  actively  engaged  in  a  research  effort 
which  will  lead  to  the  development  of  a  high-power,  Nd:YAG-based,  "pulse  burst"  laser  system. 
The  current  system  is  capable  of  creating  a  "train"  of  on  the  order  of  30-40  high  power  pulses,  at 
a  burst  frequency  of  up  to  1  MHz.  The  number  and  interpulse  spacing  of  the  individual  pulse  is 
flexible,  and  individual  pulse  energies  on  the  order  of  lO's  to  lOO's  of  mJ's  are,  ultimately, 
anticipated. 
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The  pulse  burst  laser  system  illustrated  in  Fig.  6  consists  of  the  following  four  elements: 

1.  Monolithic,  single-frequency,  diode-pumped,  cw  NdrYAG  laser. 

2.  Pulse  sheer. 

3.  Two-stage,  multi-pass  preamplifier. 

4.  Two-stage  power  amplifier. 

The  single-frequency  Nd:YAG  laser  is  a  Lightwave  Electronics  model  120-03,  with  cw  power  of 
approximately  20  mWatts.  The  output  of  this  laser  is  polarization-coupled  into  preamplifier  #1, 
which  is  a  four-pass,  10  Hz  repetition  rate,  flashlamp-pumped,  7  mm  diameter  by  110  mm  long 
Nd:YAG  rod.  The  small  signal,  single  pass  gain  has  been  measured  to  be  approximately  lOx,  so 
that  the  output  of  preamplifier  #1  is  a  single,  approximately  100  microsecond  duration,  200  Watt 
peak  power  pulse,  repeated  at  a  repetition  rat  e  of  10  Hz.  This  100  microsecond  pulse  is  formed 
into  a  "train"  using  a  pulse  sheer,  custom  built  by  Medox,  Inc.,  Ann  Arbor,  Michigan.  The  pulse 
sheer  consists  of  a  pair  of  fast  electro-optic  Pockel  cells,  and  is  capable  of  creating  an  arbitrary 
number  of  pulses,  with  minimum  individual  pulse  duration  of  10  nsec,  and  minimum  inter-pulse 
period  of  1  micro-second.  The  output  of  the  pulse  sheer  is  spatially  filtered  and  double-passed 
through  preamplifier  #2.  To  date,  we  have  obtained  10^  net  power  gain  out  of  this  system,  with 
pulse  trains  consisting  of  10-20  pulses,  and  individual  pulse  durations  between  10-20  nsec.  This 
corresponds  to  individual  pulse  energies  between  200-400  microjoules. 

The  above  described  system  has  existed  in  our  laboratory  for  approximately  six  months,  and  its 
operation  is  relahvely  straight-forward.  The  number,  spacing,  and  duration  of  the  pulses  within 
the  train  can  be  conrinuously  varied  by  flipping  a  few  thumbwheel  switches  on  the  sheer  timing 
electronics.  At  the  two-stage  power  amplifier  has  recently  become  available  and  we  are 
currently  in  the  process  of  incorporating  it  into  the  system. 


An  important  feature  of  the  laser  system  is  that  the  very  narrow  linewidth  makes  it  a  viable 
source  for  Filtered  Rayleigh  Scattering  (FRS)  measurements.  In  this  case,  the  frequency  is  tuned 
so  that  the  second  harmonic  overlaps  an  absorption  line  of  iodine  molecular  vapor.  When  an 
iodine  cell  is  placed  in  front  of  the  camera  system,  background  scattering  at  the  laser  frequency  is 
removed  from  the  image.  Light  scattered  from  the  flowing  gas,  however,  is  frequency  shifted 
due  to  the  Doppler  effect,  passes  through  the  cell,  and  is  imaged  by  the  camera.  In  this  manner, 
the  time  evolution  of  complex  flows  near  walls  an  in  ducts  can  be  imaged,  flow  velocity  and 
temperature  may  also  be  measured  if  the  cw  injection  laser  is  augmented  so  that  it  can  be  rapidly 
tuned,  leading  to  a  frequency-chirped  pulse  burst.  An  example  of  a  Filtered  Rayleigh  image  of  a 
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complex  flow  field  inside  a  duct  is  shown  in  Fig.  7.  In  this  case,  the  scattering  is  from 
condensate  ice  crystals  in  a  Mach  3  flow  entering  an  inlet  which  has  a  flat  top  and  bottom,  and 
side  fins  turned  11°  to  the  flow.  The  image  shows  individual  cross  sections  normal  to  the  flow 
taken  at  1/lOth  of  an  inch  intervals,  passing  through  a  crossing  shock  region.  In  these  images, 
boundary  layer  structure  and  shock  curvature  are  clearly  evident.  A  similar  image  taken  with  a 
pulse-burst  laser  swept  with  the  flow  could  follow  the  evolution  of  an  individual  boundary  layer 
structure  through  the  crossing  shock  region. 

The  results  to-date  on  our  pulse-burst  system  have  been  accepted  for  presentation  at  the  34th 
AIAA  Aerospace  Sciences  Meeting  in  Reiio,  Nevada  (January  1996). 
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Figure  1. 


Snosophief"^^*°"’  Comparison  of  Wind  Tunnel 
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Schematic  Diagram  of  Pulse-Burst  Laser  System. 


Cross-sectional  Filtered  Rayleigh  images  taken  at  sequential  0.1  inch 
intervals  moving  downstream  through  the  shock  crossing  region  of  an 
inlet  model. 


APPENDIX  A  —  Reprints 


R. B.  Miles,  G.L.  Brown,  W.R.  Lempert,  R.  Yetter,  GJ.  Williams,  Jr., 

S. M.  Bogdonoff,  D.  Natelson,  and  J.R.  Guest,  "Radiatively  Driven  Hypersonic 
Wind  Tunnel,"  AIAA  Journal,  Vol.  33,  No.  8,  August  1995,  pp.  1463-1470. 


S.  Macheret,  G.  Williams,  G.  Comas,  C.  Meinrenken,  W.  Lempert,  and  R.  Miles, 
Energy  Addition  and  Thermalization  Issues  in  a  Radiatively-Driven  Hypersonic 
Wind  Tunnel,"  Paper  #AIAA  95-2142,  30th  AIAA  Thermophysics  Conference, 
June  19-22, 1995,  San  Diego,  CA. 


11 


JUUKINAl. 

.  Vol.  33.  No.  8.  August  1995 
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This  paper  outlines  a  new  approach  to  the  design  of  hypersonic  wind  tunnels.  This  new  approach  is  motivated  by 
fundamental  limitations  of  conventional  iscnlropic  expansions  that  arise  from  the  very  high  temperatures  required 
to  achieve  the  necessary  enthalpy  for  hypersonic  flow.  These  high  temperatures  lead  to  excessive  throat  degradation 
and  contaminated  air  in  the  test  section.  The  consequence  is  that  the  run  times  of  conventional  facilities  must  be 
short,  and  tests  arc  conducted  in  “air”  of  unknown  composition  containing  exceedingly  high  concentrations  of 
radical  and  supcrcquilibrium  species  such  as  NO  as  well  as  ablated  plenum  and  throat  material.  The  radiatively 
driven  wind-tunnel  approach  takes  advantage  of  the  real  gas  properties  of  air  to  achieve  high  enthalpy  at  low 
temperature  in  the  plenum,  thus  minimizing  throat  degradation  and  suppressing  the  formation  of  unwanted 
species  in  the  plenum.  Additional  energy  is  radiatively  added  downstream  of  the  throat  in  the  expansion  section 
to  achieve  the  desired  test  conditions.  The  temperature  of  the  air  is  kept  low  throughout  the  expansion,  so  that 
the  formation  of  superequilibrium  species  and  radicals  is  kept  to  a  minimum.  Radiative  sources  that  couple  to  air 
include  high-power  lasers  and  microwave  devices.  A  onc-dimcnsional  model  including  optical  coupling  is  developed 
using,  as  an  example,  an  OT  laser  coupled  to  the  naturally  occurring  CO2  in  air. 


Introduction 

HE  testing  and  evaluation  of  hypersonic  airplane  components 
and  air-breathing  propulsion  systems  arc  seriously  limited  by 
the  availability  of  satisfactory  ground  test  facilities.  At  Mach  num¬ 
bers  higher  than  approximately  Mach  10,  tests  under  true  simu¬ 
lated  flight  conditions,  including  correct  atmospheric  temperature 
and  density,  cannot  be  done  with  current  or  projected  facilities.  This 
limitation.occurs  because,  in  the  conventional  iscntropic  expansion- 
type  wind-tunnel  facility,  temperatures  in  excess  of  3000  K  arc  re¬ 
quired  to  achieve  static  temperatures  on  the  order  of  200  K  or  higher 
associated  with  true  flight  conditions.  At  these  temperatures,  the  air 
cannot  be  contained  for  long  periods  of  time  due  to  materials  lim¬ 
itations,  and  significant  throat  erosion  occurs  during  wind-tunnel 
operation.  In  addition,  radical  and  superequilibrium  species  such  as 
NO  arc  formed  at  high  temperature  and  are  frozen  into  the  flow- 
field  in  the  expansion,  contaminating  the  flow.  The  solution  to  these 
problems  has  been  to  operate  with  very  short  run  times  to  pre¬ 
serve  the  throat  and  to  accept  whatever  contamination  is  present  in 
the  airflow.  This  short  time  operation  severely  limits  the  validity 
of  component  tests,  and  air  contamination  calls  into  question  any 
propulsion  measurements. 

In  this  paper  we  examine  a  new  approach  to  the  design  of  a  hyper¬ 
sonic  facility.  This  approach  leads  to  the  possibility  of  long  run  times 
with  properly  constituted  air  (no  significant  contamination)  at  static 
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temperatures  and  pressures  characteristic  of  flight  conditions.  Tlie 
basic  concept  is  to  maintain  the  air  at  moderate  temperature  but  at  ul- 
Irahigh  pressure  in  the  plenum,  generating  a  cold,  high-densify  flow 
at  the  throat.  After  an  initial  expansion  to  achieve  supersonic  flow, 
additional  energy  is  radiatively  coupled  into  the  air  as  it  expands 
through  the  nozzle,  further  accelerating  it  while  simultaneously  in¬ 
creasing  the  enthalpy  and  entropy.  A  final  expansion  to  high  Mach 
number  brings  the  flow  to  the  test  conditions.  With  this  approach, 
the  air  temperature  is  maintained  low  throughout  the  expansion  so 
that  thermally  induced  reactions  and  throat  erosion  arc  minimized. 

An  important  aspect  of  this  approach  is  that  the  ultrahigh  pressure 
air  in  the  plenum  region  has  very  high  density  and  can  no  longer 
be  treated  as  an  ideal  gas.  The  associated  real  gas  effects  cause 
the  enthalpy  to  become  pressure  dependent,  and  so  very  high  en¬ 
thalpies  can  be  achieved  at  modest  temperatures.  In  addition,  the 
speed  of  sound  increases  significantly,  reaching  values  more  typical 
of  the  speed  of  sound  in  liquid  than  in  gas-phase  air.  As  a  conse¬ 
quence,  the  kinetic  energy  of  the  fluid  passing  through  the  throat  at 
Mach  1  is  much  higher  than  would  be  possible  for  an  ideal  gas  at  a 
similar  temperature. 

A  feature  of  this  approach  is  that  many  parameters  may  be  varied 
to  achieve  the  desired  flow.  For  example,  the  absorption  of  the  air 
depends  on  its  density,  temperature,  and  pressure,  as  well  as  on  the 
wavelength  of  the  particular  radiative  source  chosen.  A  change  in 
the  absorption  will  affect  the  length  of  the  tunnel,  which  directly 
impacts  the  time  the  air  spends  at  high  temperature  and,  thus,  the 
chemical  dynamics,  and  so  the  choice  of  the  wavelength  can  be  op¬ 
timized  for  the  minimization  of  oxides  of  nitrogen.  Other  variables 
include  the  plenum  pressure,  the  plenum  temperature,  the  expan¬ 
sion  profile,  and  the  energy  addition  profile.  More  flexibility  may 
be  achieved  by  using  multiple  radiation  sources,  each  optimized  for 
a  particular  portion  of  the  flow.  A  practical  limit  for  the  run  time 
will  be  determined  by  the  volume  of  high-pressure  air  that  can  be 
contained  and  the  total  stored  energy  of  the  radiative  source.  Run 
times  in  excess  of  0.1  s  are  the  goal  of  this  approach. 

Technical  Approach 

Overview  T  ■ 

A  convenient  way  of  viewing  the  operation  of  such  a  facility  is 
with  the  help  of  a  Mollicr  diagram  (Fig.  1).  In  this  diagram,  the 
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Fig.  1  Concept  definition:  thermodynamic  comparison  of  wind-tunnel  philosophies. 


Fig.  2  Equilibrium  mole  fractions  of  NO,  NO2,  and  N2O  in  air  at  1 
and  50,000  atm.  For  1  atm,  the  ideal  gas  equation  of  state  was  applied, 
and  for  50  atm,  the  Becker-  Kistiakowsky-Wilson  equation  of  state  was 
assumed.^ 

enthalpy  is  plotted  vs  the  entropy.  For  a  conventional  wind-tunnel 
facility,  the  expansion  is  assumed  to  be  isentropic,  and  so  the  static 
temperature  and  pressure  and  the  Mach  number  in  the  test  cham¬ 
ber  mandate  a  unique  temperature  and  pressure  for  the  gas  in  the 
plenum.  For  example,  for  a  Mach  16  flow  with  a  static  tempera¬ 
ture  of  200  K  and  a  static  pressure  of  10'^  atm  (point  C  on  the 
diagram),  the  plenum  must  be  held  at  a  temperature  of  6300  K  and 
at  a  pressure  of  4000  atm  (point  A).  Clearly,  at  this  temperature, 
air  is  highly  dissociated,  and  severe  problems  with  containment 
and  throat  erosion  occur.  In  contrast,  the  radiativciy  driven  wind 
tunnel  begins  with  a  plenum  pressure  on  the  order  of  30,000  atm 
but  a  plenum  temperature  on  the  order  of  1000  K.  The  point  on 
the  Mollier  diagram  (point  B)  associated  with  this  initial  condi¬ 
tion  corresponds  to  significantly  lower  entropy  and  lower  enthalpy 
than  that  of  the  isentropic  wind  tunnel.  That  remaining  enthalpy 
and  entropy  arc  then  added  downstream  of  the  throat  using  the 
radiative  source. 

Chemistry  is  an  important  consideration  in  the  design  of  the  ra¬ 
diativciy  driven  hypersonic  facility  concept.  In  the  ideal  ease,  the 
temperatures  will  be  kept  low  enough  so  that  the  chemical  kinet¬ 
ics  can  be  ignored.  In  reality,  however,  there  arc  two  regimes  in 
which  chemical  kinetics  must  be  taken  into  account.  The  first  is 
in  the  ultra-high-pressure  plenum,  where  even  at  modest  tempera¬ 
tures  the  formation  of  polyatomic  molecules  that  have  lower  specific 
volume  is  favored.  This  effect  is  illustrated  in  Fig.  2.  If  equilib¬ 
rium  codes  arc  extrapolated  to  50,000  atm  at  1000  K,  the  predicted 
equilibrium  mole  fraction  of  NO2  increases  by  two  orders  of  mag¬ 
nitude  above  its  one  atmosphere  equilibrium  concentration.  The 
mole  fraction  of  NO,  on  the  other  hand,  is  unaffected  since  it  is  a 


diatomic  molecule  and  occupies  the  same  molar  volume  as  oxygen 
and  nitrogen,  which  arc  both  diatomic.  At  temperatures  on  the  order 
of  1500  K  or  below,  the  rate  of  formation  of  those  species  is  very 
slow,  even  at  high  pressure,  and  so  a  facility  pulsed  for  times  on  the 
order  of  a  second  or  less  may  never  reach  chemical  equilibrium  in 
the  plenum. 

The  second  regime  of  concern  is  in  the  region  of  heat  addition 
where  the  temperatures  increase  substantially,  and  the  formation 
rate  of  superequilibrium  species  is  nonncgligiblc.  Since  the  flow 
is  moving  at  high  speed  (several  kilometcrs/sccond).  the  residence 
time  of  the  molecules  in  the  high  temperature  regime  is  relatively 
short  (milliseconds),  and  so  the  kinetic  rate  of  formation  of  various 
molecular  species  becomes  an  important  parameter.  In  particular,  it 
is  important  to  realize  that  the  slow  kinetic  rate  of  formation  of  the 
nitric  oxide  molecule  may  suppress  the  nitric  oxide  concentrations 
for  short  times  even  though  the  temperature  may  be  high. 

One  can  generate  an  estimate  of  the  temperature  at  which  the  ad¬ 
ditional  enthalpy  must  be  added  downstream  by  recognizing  that  the 
total  enthalpy  (kinetic  energy  plus  static  enthalpy)  and  entropy  arc 
determined  by  the  desired  test  conditions.  Therefore,  we  can  iden¬ 
tify  an  enthalpy  deficit  A/z,  which  is  the  difference  between  the  total 
enthalpy  in  the  test  chamber  and  the  enthalpy  in  the  plenum,  and 
an  entropy  deficit  Aj,  which  is  the  difference  between  the  entropy 
of  the  gas  in  the  lest  chamber  and  the  entropy  in  the  plenum.  In  the 
conventional  isentropic  facility,  the  enthalpy  and  entropy  deficits 
arc  zero.  Assuming  an  isentropic  expansion,  the  total  enthalpy  in 
the  test  chamber  is  equal  to  the  plenum  enthalpy  in  the  conventional 
facility,  and  so  the  enthalpy  deficit  is  just  the  enthalpy  difference 
between  the  plenum  enthalpy  of  the  conventional  wind  tunnel  and 
the  plenum  enthalpy  of  the  radiatively  driven  wind  tunnel,  as  shown 
in  Fig.  1.  Similarly,  the  entropy  deficit  is  the  difference  in  plenum 
entropies  as  shown  in  the  figure.  When  energy  is  added  downstream 
of  the  throat,  the  ratio  of  these  two  deficits  (Afi/^s)  gives  an  “ef¬ 
fective”  temperature  T^,  This  is  the  effective  temperature  at  which 
the  additional  enthalpy  must  be  introduced  downstream  of  the  throat 
to  achieve  the  proper  entropy. 

To  minimize  the  formation  of  oxides  of  nitrogen  and  radicals  in 
the  nozzle  section,  it  is  desirable  to  keep  this  effective  tempera¬ 
ture  as  low  as  possible.  This  is  accomplished  by  either  decreasing 
Ah  (increasing  the  enthalpy  in  the  plenum)  or  increasing  As  (de¬ 
creasing  the  entropy  in  the  plenum).  From  Fig.  I  it  is  apparent  that 
increasing  the  enthalpy  in  the  plenum  can  be  achieved  by  increasing 
the  temperature  or  by  increasing  the  pressure.  Decreasing  the  en¬ 
tropy  is  accomplished  by  increasing  the  pressure  or  by  decreasing 
the  temperature.  It  is  evident  that  high  pressure  is  needed  to  keep 
the  effective  temperature  low.  Figure  3  shows  the  effective  tem¬ 
perature  as  a  function  of  plenum  pressure  for  various  plenum 
temperatures  needed  to  achieve  a  Mach  16  flow  at  200  K 
and  10"^  atm. 
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p  and  r  are  defined  in  terms  of  the  critical  values  for  the  substance  to  be  modeled: 
l<TcP  ,  Tc 


RTcP  Tc 

r  =  and  r  =  — 

Pc  r 

where  the  critical  parameters  for  air  arc  approximated  as 


Tc  =  132.5  K 


Pc  =  3.77  MPa 


Pc  =  343.3  kg/nr 
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Fig.  3  Effective  temperature  vs  plenum  pressure  for  various  plenum 
temperatures. 


Balanced  against  this  effective  temperature  are  the  .  maxi¬ 
mum  temperature  reached  in  the  expansion  and  the  temperature 
in  the  plenum.  In  general,  it  is  these  temperatures  and  the  lime  over 
which  the  flow  is  at  these  temperatures  that  determine  the  concentra¬ 
tion  of  nitric  oxide  in  the  test  region.  Since  the  effective  temperature 
is  an  average  temperature  at  which  the  energy  is  added,  if  a  signif¬ 
icant  amount  of  energy  is  added  below  the  effective  temperature, 
then  an  approximately  equal  amount  must  be  added  at  a  tempera¬ 
ture  higher  than  the  effective  temperature  to  achieve  the  proper  test 
conditions.  The  maximum  heating  rate  downstream  of  the  throat  is 
at  constant  pressure,  since  an  increase  in  pressure  will  cause  the 
flow  to  choke.  Thus,  to  minimize  the  peak  temperature,  the  flow 
should  be  heated  isobarically  until  it  reaches  a  temperature  some¬ 
what  higher  than  the  effective  temperature,  and  then  the  remainder 
of  the  energy  added.  In  the  models  presented  here,  we  have  cho¬ 
sen  an  initial  expansion  to  Mach  2  before  energy  is  added  to  avoid 
instabilities  near  Mach  1.  Because  of  real  gas  effects,  in  the  initial 
constant  pressure  heating  region  the  Mach  number  does  not  drop 
significantly,  even  though  the  temperature  increases  by  a  factor  of 
3  or  more. 


Flow  Model 

A  more  accurate  description  of  the  operation  of  the  wind  tun¬ 
nel  can  be  developed  using  a  one-dimensional  real  gas  model.Thc 
governing  equations  are  the  continuity,  momentum,  and  energy 


equations  that  arc  derived  from  a  control  volume  analysis  with  heat 
input  q,  expressed  in  joules  per  unit  mass  per  unit  length, 

6A  du  dp 

pu- — h  p  —  A  +  — =  0  (1) 

cLc  djc  dx 


dP\  dp  fdP\  dT  du 

dp )^±K  \(ir  ~ 


f  dh\  (ip  f  dh\  dT  (iH  q 

jj-dT d7 

where  p  is  tlie  density,  h  is  the  specific  enthalpy,  and  //  is  llie  velocity. 
The  equation  of  state  is^ 

P  =  Zip,  T)pRT  (4) 

where  the  compressibility  factor  is  given  as 

10  13 

Z  =  1  +  r/)(r)  +  /-^  ^  A'(/)t''"" 

1=7  1  =  11 

+  r''A:(14)r  +  r^[A:(15)T^  +  A'(16)t^] 

xr‘A'(17)r^  +  r’[/:(18)T^  +  K09)x^\ 

+  r'^K(20)r^ 

-  r^[K  +  K(,22)x'^] 

+  r^[Ki23)x^  +  K(24)x^] 

+  r''lK(25)x^  +  K(.26)xU 

+  r®[/:(27)t^  + A’(28)t5]  ^ 

+  r“>[A:(29)T’  + A:(30)t‘'] 

_  +r‘^[/(:(31)T^  + A’(32)t'’ +  /:(33)r*]_ 

with  r  being  a  dimensionless  density  based  on  ilic  critical  tempera¬ 
ture  Tc  and  the  critical  pressure  Pc, 

RTc 

r-P-  (6) 

and  T  being  a  dimensionless  temperature, 

T. 

r  =  ^  (7) 

The  constants  K{i)  arc  listed  in  Table  1.  To  find  the  other  thermo¬ 
dynamic  state  variables,  the  ideal  gas  specific  heat,  CyCT),  is  used  to 
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find  the  low-density  (p  =  p^)  (ideal  gas)  conditions  at  temperature 
7.  and  then  the  real  gas  equation  is  used  to  integrate  at  constant  T 
to  the  desired  density.^  Therefore,  for  enthalpy  and  entropy, 


where  =  pf.{RTc/ Pc)  is  effectively  zero  for  the  lower  limits  of 
the  integrals,  and  /z”  contains  the  low-density  specific  enthalpies  of 
oxygen  and  nitrogen  and  contains  the  low-density  specific  entropies 
of  oxygen  and  nitrogen  plus  the  entropy  of  mixing. 

This  model  has  been  found  to  give  an  accurate  reproduction  of 
the  Arnold  Engineering  and  Development  Center  (AEDC)  Mollier 
diagram  that  is  derived  from  tabulated  data,  and  it  is  consistent 
with  similar  real  gas  models  published  by  the  Soviet  Bureau  of 
Standards.^~^  This  equation  of  state  is  solved  using  a  finite  differ¬ 
ence  solution  to  generate  an  aerodynamic  model  to  describe  the  flow 
evolution  with  heating  in  the  nozzle  section.  The  model  does  not  ac¬ 
count  for  dissociation  or  high  temperature  and  has  been  developed 
based  on  nitrogen  data  only  up  to  10,000  atm  pressure  at  tempera¬ 
tures  between  50  and  1200  K.  Extrapolations  to  higher  pressure  and 
temperature  arc  smooth  but  have  not  been  experimentally  verified. 
The  application  of  this  model  to  air  rather  than  nitrogen  involves 
a  single  substance  approximation  for  the  critical  pressure  and  tem¬ 
perature.  To  our  knowledge  there  are  no  experimental  data  for  air 
above  4550  atm.^ 

This  model  was  tested  in  several  manners.  First,  with  no  energy 
addition  and  at  constant  area,  calculated  values  of  the  density,  tem¬ 
perature,  and  velocity  were  checked  to  confirm  that  these  quantities 
remained  constant  throughout  the  flow.  A  second  check  was  done 
modeling  a  diverging  conical  profile  with  initial  conditions  corre¬ 
sponding  to  the  Mach  number  just  greater  than  1  and  no  energy 
addition.  Calculated  p/p\  T/T\  P/P\  and  A/A*  profiles  pro¬ 
duced  by  the  model  were  compared  with  tabulated  solutions  for  an 
iscntropic  expansion.  There  was  no  perceptible  deviation  from  the 
isentropic  case.  A  third  verification  used  the  exactly  solvable  prob¬ 
lem  of  energy  addition  into  a  flow  with  a  constant  cross-sectional 
area.  The  exact  solution  in  this  case  is  given  by  the  expression 


d/z  = 


u^a 


(Is 


where  or  is  the  isobaric  expansion  coefficient, 


(10) 


1 


Of  =  — 

u 


(H) 


and  V  is  the  specific  volume  (1/p)  .  Both  a  and  the  specific  heat 
Cp  were  calculated  numerically,  and  the  ratio  ((lli/6x)/ids/dx)  was 
normalized  by  Eq.  (10)  and  checked  to  be  sure  that  the  value  re¬ 
mained  unity  as  the  flow  evolved  down  the  constant  area  duct.  The 
result  was  accurate  to  one  part  in  10\  independent  of  p,  7.  u,  and 
q  for  appropriately  small  step  sizes.  The  computation  uses  as  an  in¬ 
dependent  variable  cither  the  area  profile  A(jc),  the  pressure  profile 
p(x),  or  the  temperature  profile  T(x),  TIic  energy  deposition  rate 
q(x)  can  cither  be  specified,  or  it  can  be  derived  from  the  absorp¬ 
tion  characteristics  of  the  air,  which  arc  modeled  as  a  function  of 
temperature  and  pressure. 


Energy  Addition 

The  performance  of  the  wind  tunnel  depends  on  both  the  plenum 
conditions  and  on  the  mechanism  for  energy  addition.  Using  current 
technology,  the  maximum  plenum  pressure  achievable  is  estimated 
to  be  between  20,000  and  30,000  atm.  Energy  must  be  added  to 
the  moving  flow,  and  so  radiative  coupling  is  an  obvious  approach, 
^is  coupling  must  be  accomplished  via  absorptive  energy  transi¬ 
tions  that  occur  in  air.  The  very  high  pressure  of  the  air  just  down¬ 
stream  of  the  throat  will  lead  to  significant  pressure  dependence 


Conventional  Hypersonic  Tunnel 


Radiation-driven  Hypersonic  Tunnel 


Ring  of  Laser  Beams 


Fig.  4  Schematic  comparison  of  wind-tunnel  philosophies. 


(both  pressure  narrowing  and  pressure  broadening)  of  absorption 
features,  and  at  these  high  pressures,  pressure-induced  absorption 
phenomena  are  often  dominant.  Transitions  occurring  in  the  infrared 
region  arc  accessible  to  various  high-power  lasers,  and  those  in  the 
microwave  region  can  be  pumped  by  high-power  gyrotron  devices. 
Of  particular  interest  are  pressure-induced  transitions  in  oxygen  in 
the  near  infrared  that  overlap  diode  lasers  and  diode-pumped,  solid- 
state  lasers;  CO2  transitions  that  overlap  hydrogen  fluoride  (HF)  and 
reactor  pumped  xenon  laser  lines;  H2O  transitions  that  overlap  the 
chemical  oxygen  iodine  laser  (COIL);  and  the  oxygen  microwave 
transition  at  60  GHz.  Other  transitions  to  be  explored  include  H2O 
and  CO2  lines  that  overlap  CO2  and  CO  lasers  and  pressure-induced 
oxygen  and  nitrogen  transitions  throughout  the  infrared.  Ultravio¬ 
let  transitions,  particularly  those  in  oxygen,  might  also  be  explored 
with  high-power  light  sources  and,  possibly,  with  ultraviolet  lasers. 

As  a  test  case,  the  possibility  of  adding  energy  optically  with 
high-power  HF  lasers  was  studied.  Although  the  calculation  is  one 
dimensional,  a  possible  tunnel  geometry  is  sketched  in  Fig.  4.  Tlie 
lasers  are  assumed  to  enter  the  tunnel  from  downstream,  passing 
through  the  test  section  and  into  the  nozzle.  Tlicy  arc  arranged  in 
a  ring  so  as  not  to  interfere  with  the  model.  Since  most  of  the  ab¬ 
sorption  occurs  where  the  air  is  at  high  pressure,  this  configuration 
has  the  potential  of  efficiently  coupling  energy  into  the  flow  just 
downstream  of  the  throat.  To  compute  the  flow,  an  initial  guess 
of  the  energy  deposition  parameter  q  through  the  nozzle  is  given, 
and  the  model  is  integrated  in  the  forward  direction  to  solve  for 
the  air  density  and  temperature  through  the  expansion,  and  then 
in  the  backward  direction  (from  the  test  section  back  toward  the 
throat)  to  give  a  second  iteration  on  the  energy  addition  profile 
based  on  the  energy  absorbed  by  the  air  at  the  computed  temper¬ 
ature  and  density.  This  procedure  is  carried  out  repeatedly  until 
convergence  is  achieved.  Following  convergence,  the  energy  addi¬ 
tion  rate  profile  is  consistent  with  the  absorptivity  profile  of  the 
expanding  airstream. 

Although  the  thermodynamic  properties  of  the  gas  appear  to  be 
good  approximations  to  the  actual  thermodynamic  state  of  air  over 
these  ranges  of  pressure  and  temperature,  the  absorption  constant  is 
only  a  crude  estimate.  In  the  example  worked  here,  the  high-power 
HF  lasers  arc  assumed  to  be  operating  on  a  single  transition  near 
2.76  /z  (3623  cm~‘).  Absorption  occurs  through  the  small  amount 
of  CO2  that  is  present  in  air  and  reflects  the  upper  atmospheric  mole 
fraction  ratio  of  330  parts/million.  The  absorption  is  directly  into 
a  bcnding/asymmctric  stretch  (021)  combination  band,  and  absorp¬ 
tion  data  arc  taken  from  the  HITRAN  database,  extrapolated  to  high 
pressure  using  conventional  pressure  broadening,’  Since  pressure 
narrowing  effects  arc  ignored,  the  values  that  arc  used  only  serve  to 
illustrate  this  approach  to  energy  addition  and  cannot  be  interpreted 
quantitatively.  Figure  5  shows  the  absorption  profile  in  the  vicin¬ 
ity  of  the  HF  line  at  a  pressure  of  1000  atm.  To  accommodate  our 
uncertainty  in  the  absorption  line  shape  and  magnitude,  a  pressure- 
dependent  absorption  constant  equivalent  to  one-half  of  that  shown 
on  the  figure  for  the  Pi{A)  line  was  used  for  the  calculation. 
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Tabic  2  Computed  gases 


Stagnation  conditions 

Cases  1  and  2 

Temperature,  K 

1500 

Pressure,  atm 

2.51  X  10' 

Density,  kg/m^ 

1.12  X  103 

Enthalpy,  J/kg 

3.54  X  10^ 

Entropy.  J/kg  -  K 

5.56  X  10^ 

Flow  rate,  kg/s 

79.0 

Throat  conditions  (Mach  1) 

Velocity,  m/s 

1.86  X  103 

Temperature,  K 

1080 

Density,  kg/m^ 

8.17  X  103 

Area,  m^ 

5.19  X  10-3 

Pressure,  atm 

8.33  X  103 

Exit  conditions 

Case  1 

Case  2 

Mach  no. 

14.5 

15.0 

Velocity,  m/s 

4.05  X  10^ 

4.25  X  103 

Temperature,  K 

194 

201 

Density,  kg/m^ 

1.82  X  10-3 

1.76  X  10-3 

Area, 

10.7 

10.6 

Pressure,  atm 

1.00  X  10-3 

1.00  X  10-3 

Sialic  enthalpy,  J/kg 

1.90  X  103 

1.90  X  10* 

Total  enthalpy,  J/kg 

8.40  X  10* 

9.22  X  10* 

Entropy,  J/kg  •  K 

8.40  X  103 

8.44  X  103 

Equivalent  isentropic  plenum 

conditions  (from  AEDC  Mollier  diagram) 

Temperature,  K 

5550 

5800 

Pressure,  atm 

2.24  X  10^ 

2.82  X  103 

IIFUm(  Unti 

too 

90 
80 
70 
60 

Absorption  .. 

(%/m) 

40 

30 
20 
10 
0 

3500  3550  3600  3650  3700  3750  3800 

Frequency  (1/cm) 

Fig.  5  COz' absorption  (percent  per  meter)  for  air  at  1000  atm  and 
1000  K,  Vertical  lines  correspond  to  IIF  laser  lines. 


Fig.  6  MoIIicr  diagram  for  eases  I  (2000  K— solid  line,  solid  dot)  and 
2  (2500  K — dotted  line,  open  dot). 


Figure  6  shows  the  trajectories  of  eases  1  and  2  on  the  MoIIicr 
diagram.  The  points  at  the  top  of  the  diagram  show  the  initial  condi¬ 
tions  for  an  iscntropic  expansion  to  achieve  the  same  test  conditions. 
These  numbers  arc  taken  from  the  AEDC  MoIIicr  diagram  since  our 
model  is  not  designed  for  exceedingly  high  temperature. 

Figures  7a-7f  show  the  pressure,  temperature,  tunnel  radius, 
Mach  number.  How  velocity,  and  heating  profile  as  a  function  of 
the  axial  distance  downstream  for  these  two  eases.  It  is  particularly 
striking  to  note  the  large  change  in  tunnel  radius  as  a  function  of 
axial  distance.  Figures  8a  and  8b  show  an  expanded  view  of  the  ra¬ 
dius  and  the  heating  profile  over  the  first  few  meters.  Because  of  the 
very  high  density  of  the  fiuid  at  the  throat,  the  radius  at  that  location 
is  less  than  5  mm.  Apparent  discontinuities  in  these  curves  arc  due 
to  discontinuous  changes  in  the  modeled  expansion  profiles  or  en¬ 
ergy  addition  constraints.  As  the  fiow  proceeds  downstream,  energy 
addition  occurs  primarily  over  the  first  few  meters,  after  which  the 
secondary  expansion  increases  the  radius  to  on  the  order  of  1.7  m. 
The  fiuid  fiow  characteristics  at  the  throat  are  shown  in  Table  2.  Note 
that  the  velocity  at  the  throat  is  1,860  m/s.  This  sonic  speed  would 
correspond  to  an  “ideal  gas”  at  a  temperature  in  excess  of  8000  K, 
whereas  the  fiow  temperature  is  actually  1080  K, 

In  ease  2,  significantly  more  energy  is  coupled  into  the  flow  since 
it  is  heated  to  2500  K  rather  than  2000  K  before  the  isothermal 
energy  addition  begins.  This  is  rcfiected  both  in  the  test  conditions 
and  in  the  effective  temperatures.  In  case  I ,  the  effective  temperature 
is  1740  K,  the  fiow  velocity  in  the  test  section  is  4.05  km/s,  and  the 
static  temperature  is  194  K.  For  ease  2,  the  effective  temperature  is 
1970  K,  the  fiow  velocity  is  4.25  km/s,  and  the  temperature  in  the 
test  section  is  201  K. 


Two  examples  are  given  to  demonstrate  the  effect  of  different  ex¬ 
pansions.  Both  eases  begin  with  the  same  plenum  pressure,  25,120 
atm.  and  temperature,  1500  K,  and  arc  expanded  to  roughly  200  K 
and  10“^  atm.  Mach  numbers  of  14.5  and  15.0  arc  achieved  by 
eases  1  and  2,  respectively.  The  stagnation  conditions,  throat  con¬ 
ditions,  and  conditions  in  the  test  section  for  the  computed  eases 
arc  shown  in  Table  2.  The  constraint  applied  for  the  calculation  was 
initial  energy  addition  at  constant  pressure,  followed  by  isothermal 
energy  addition,  and  a  final  expansion  to  the  test  cell  that  has  ap¬ 
proximately  a  3-m  diameter.  In  the  first  case,  the  fiow  is  heated  at 
constant  pressure  to  2000  K  and  then  heated  at  constant  temperature. 
Since  the  code  is  one  dimensional,  the  length  and  shape  of  the  final 
expansion  arc  arbitrary.  The  shape  of  the  final  expansion  is  chosen 
here  to  produce  a  smooth  pressure  drop  to  10”^  atm  pressure,  and 
the  length  of  the  tunnel  is  chosen  to  be  40  m.  In  the  second  ease,  the 
flow  is  healed  at  constant  pressure  up  to  2500  K  and  then  heated  at 
constant  temperature  until  the  proper  entropy  is  reached.  The  initial 
conditions  for  the  energy  addition  portion  of  the  computation  begin 
at  Mach  2  in  both  eases,  and  so  the  flow  passes  through  the  throat 
and  expands  to  Mach  2  before  heat  addition  begins.  In  both  eases, 
the  optical  energy  is  fully  absorbed  downstream  of  the  Mach  2  po¬ 
sition  by  allowing  sufficient  separation  between  the  nozzle  location 
and  the  energy  addition  regime,  and  so  there  is  no  truncation  error 
associated  with  these  initial  conditions. 


Flow  Chemistry 

An  important  design  goal  for  this  wind-tunnel  facility  is  to  match 
the  molecular  concentrations  that  arc  found  at  the  flight  altitudes 
simulated.  Particular  difficulty  arises  with  regard  to  the  molecu¬ 
lar  concentration  of  nitric  oxide,  since  that  species  is  formed  in 
abundance  at  high  temperature  and  tends  to  freeze  into  the  fiow. 
Nitric  oxide  concentrations  in  the  upper  atmosphere  range  from  on 
the  order  of  50  parts/billion  up  to  80  km  (260,000  ft)  and  then  in¬ 
crease  to  a  maximum  of  approximately  200  parts/million  at  1 10 
km  (360,000  ft).®  These  mole  fractions  correspond  to  equilibrium 
temperatures  ranging  from  approximately  600  to  1200  K,  according 
to  the  equilibrium  diagram  shown  in  Fig.  2.  An  important  feature 
of  this  facility  is  that,  even  though  the  air  must  reach  temperatures 
higher  than  these  numbers,  the  time  at  which  it  remains  at  those  tem¬ 
peratures  is  short,  and  therefore  the  nitric  oxide  concentrations  arc 
significantly  suppressed  below  their  equilibrium  levels.  This  occurs 
both  in  the  plenum  where,  even  though  the  pressure  is  very  high, 
the  temperature  is  low,  and  so  the  kinetic  formation  rale  of  nitric 
oxide  is  on  the  order  of  seconds,  and  in  the  expansion,  where  the 
temperatures  are  higher  but  the  limes  for  formation  are  only  on  the 
order  of  milliseconds. 

■To  compute  the  NO  concentrations,  we  have  applied  an  elemen¬ 
tary  chemical  kinetic  model  based  on  23  principal  air  reactions  and 
the  real  gas  Beckcr-Kisliakowsky- Wilson  equation  of  state.  These 
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Fig.  7  T^.nncl  parameters  for  cases  1  (2000  K-soIid  line)  and  2  (2500  K-dotted  line). 


reactions  and  the  pressure-dependent  rates  are  given  in  Table  3 

concentrations  for  the  two 
no!e^ haHi"''  TK 9  and  1 0.  It  is  important  to 

K  IS  approximately  lO'^  (}%).  Here  we  see  that  in  case  1  the  for- 
dicmicaff"""^^'^*'  dramatically  suppressed  by  the  slow 

clcentration  Ll  f  ^  'ong  Period  of  time,  and  the  nitric  oxide 
concentration  comes  close  to  its  equilibrium  value  These  curves 

Sfa  scenarios  for  energy  addition  will 

a  significant  impact  on  the  level  of  nitric  oxide  formed  The 

bSloirneX" 


In  this  model,  the  time  delay  between  the  excitation  of  the  CO, 
molecule  and  the  deposition  of  energy  into  the  kinetic  motion  of 
the  now  is  modeled  assuming  a  single  vibrational  temperature 
Although  the  exact  mechanisms  for  this  energy  transfer  remain 
0  be  explored,  it  is  expected  that  the  COj  molecule  will  rapidly 
transfer  its  energy  to  nitrogen  by  exciting  the  vibrational  mode  of 
he  nitrogen  molecule.  At  low  pressure,  vibrationally  excited  ni¬ 
trogen  has  a  very  long  relaxation  time.  At  high  pressures,  how¬ 
ever,  nitrogen  rapidly  transfers  its  excitation  to  oxygen,  which  is 
rapidly  thermalizcd.  Accurate  modeling  of  both  the  collisional  en¬ 
ergy  transfer  from  CO2  to  nitrogen  and  the  nitrogen  relaxation  mech¬ 
anism  wfil  be  important  since  many  photons  must  be  transferred 
to  each  CO2  molecule  to  achiev:  ;hc  required  heating  rates.  Tlic 


MILLS  fiTAL 


J469 


Tabic  3  High-tcmpcraturc  air  reaclion  mechanism^ 


Reaction^ 

A 

n 

E 

fH) 

^'0 

a 

b 

1 

OjC+M)  =  20(+M) 

0.15  X  10‘2 

0.0 

118,000.0 

0.98  X  10“ 

-2.5 

118.000.0 

2 

0  +  02(+M)  =  03(+M) 

0.17  X  10‘^ 

0.0 

- 0.0 

0.18  X  10®® 

-2.8 

_ f)  () 

0.65 

3 

0  -1-  O3  =  2O2 

0.12  X  lO'"* 

0.0 

- 4.5 

_ 

_ 

_ 

4 

CO  +  0  +  M  =  CO2  +  M 

— 

— 

— 

0.25  X  10‘“ 

0.0 

-4,541.0 

5 

CO+O2  =  CO2  -f  0 

0.25  X  10'^ 

0.0 

47,690.0 

_ 

_ 

6 

NO(+M)  =  N  +  0(+M) 

0.15  X  10‘^ 

0.0 

148,000.0 

0.15  X  10“ 

0.0 

148.400.0 

7  . 

NO  +  0(+M)=N02(+M) 

0.13  X  10'* 

-0.75 

0.0 

0.47  X  10®® 

-2.87 

0 

0 

0.95 

1.0  X  10“^ 

8 

NO  +  0  =  O2  +  N 

0.18  X  10'® 

1.0 

38,760.0 

_ 

.... _ 

_ _ 

9 

2NO  =  N2  +  O2 

0.13  X  10'^ 

0.0 

75,630.0 

_ 

_ 

10 

NO2  +  0  =  O2  +  NO 

0.39  X  10'^ 

0.0 

-238.4 

_ 

_ 

11 

NO2  +  0(+M)  =  N03(+M) 

0.13  X  10'“ 

0.0 

0.0 

0,15  X  10®^ 

-4.08 

1.467.0 

0.79 

2.5  X  10"^ 

12 

NO2  +  NO  =  N2O  +  Q, 

0.10  X  10'^ 

0.0 

60,000.0 

— 

_ 

_ 

_ 

_ 

13 

2NO2  =  NO3  +  NO 

0.96  X  10'® 

0.73 

20,920.0 

_ 

_ 

_ 

14 

2NO2  =  2NO  +  O2 

0.16  X  10‘® 

0.0 

26,120.0 

_ 

_ 

_ 

15 

N20(+M)  =  N2  0(+M) 

0.13  X  10'^ 

0.0 

59,610.0 

0.72  X  10'® 

-0.73 

62,790.0 

1.17 

1.25  X  10"^ 

16 

N2O  -i-  0  =  N2  +  O2 

0.10  X  10'5 

0.0 

28.020.0 

_ 

_ 

_ 

17 

N2O  4*  0  =  2NO 

0.66  X  10'“ 

0.0 

26,630.0 

_ 

_ 

_ 

18 

N2O  -1-  NO  =  N2  +  NO2 

0.10  X  10'^ 

0.0 

49.675.0 

_ 

_ 

_ 

19 

N  NO  =  N2  +  0  • 

0.33  X  10'® 

0.3 

0.0 

_ 

_ 

_ 

_ 

20 

N  N2O  =  N2  4-  NO 

0.10  X  10'“ 

0.0 

19,870.0 

_ 

_ _ _ 

_ 

_ 

_ 

21 

N  4*  NO2  =  2NO 

0.40  X  10'® 

0.0 

0.0 

_ 

_ 

_ 

_ 

_ 

22 

N  4-  NO2  =  N2O  4-  0 

0.50  X  10'® 

0.0 

0.0 

_ 

_ 

_ 

_ 

_ 

23 

N4-NO2  =  N2  4-  O2 

0.10  X  10'® 

0.0 

0.0 

— 

— 

— 

— 

— 

“Unils  arc  cm^-inolc-s-cal,  k=z  AX'*  C7^p{~-Ea| RT). 

*^ror  prcssurc-dcpcndcnl  rcaciions  with  the  total  concentration.  M  listed  in  parentheses,  fall-off  behavior  is  included  as  follows: 


c- 

Ao  +  (W[«])  1  ^ 


ll  +  (loE(to  x(iW]/*„)l2) 
*0  =  AoT"o  cxp(-£o//J7-) 

=  A7''’cxp(-£//?r) 


Fc  =  a  ~  bT 


a)  Radius 


Axial  Position 
(m) 


b)  Heating  profile 


Axial  Position 
(m) 


Fig.  8  Expanded  scale  wind-tunnel  parameters  for  eases  1  (2000  K— solid  line)  and  2  (2500  K— doited  line). 


scale  of  this  energy  deposition  problem  can  be  better  appreciated 
by  realizing  that  the  enthalpy  deficits  in  eases  1  and  2  arc  4.9  and  ' 
5.7  million  J/kg,  respectively.  This  corresponds  to  approximately 
10,000  photons  per  CO2  . molecule,  all  of  which  must  be  absorbed 
in  a  few  milliseconds,  and  so  collisional  transfer  rates  shorter  than 
0. 1  /xs  arc  required.  In  a  large-scale  facility  such  as  this,  a  flow  rate 
on  the  order  of  80  kg/s  will  be  necessary.  Thus,  a  total  energy  de¬ 
position  rate  of  400-450  MW  must  be  achieved  over  the  run  time 
of  the  tunnel.  This  involves  a  significant  scaling  of  candidate  laser 
sources  and  the  use  of  multiple  lasers.  Tlinncl  operation  for  0.1  s 
will  require  a  total  energy  of  50  MJ  to  be  deposited  into  the  flow.  By 
comparison,  this  is  approximately  the  amount  of  energy  stored  in 


20  gal  of  gasoline  and  is  well  within  the  current  total  energy  capa¬ 
bility  of  existing  high-power  hydrogen  fluoridc/deutcrium  fluoride 
(HF/DF)  lasers. 

Alternate  Approaches 

The  use  of  an  HF  laser  to  achieve  energy  addition  via  a  CO2  transi¬ 
tion  has  the  advantage  that  very  high-power  HF  and  DF  lasers  have 
already  been  developed  through  Department  of  Defense  support. 
For  example,  the  MIRACL  laser  currently  is  capable  of  operating 
at  power  levels  in  excess  of  I  x  10^  W.  Further  scaling  of  that 
laser  to  power  levels  on  the  order  of  300  x  10^  W  will  constitute 
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Axial  distance 
(m) 

Fig.  9  Nonequilibrium  concentrations  for  ease  1. 


Axial  distance 
(m) 


Fig.  10  Nonequilibrium  concentrations  for  case  2. 


a  major  engineering  challenge.  Other  candidate  laser  sources  for 
energy  transfer  include  the  nuclear  reactor  pumped  Xc  laser  oper¬ 
ating  at  2.76  fi,  the  free  electron  laser  with  its  variable  frequency 
operation,  the  chemical  oxygen  iodine  laser  that  runs  at  1.4  /x,  high 
efficiency  diode  lasers  operating  in  the  near  infrared,  the  carbon 
dioxide  laser  in  the  10-;x  region,  and  the  carbon  monoxide  laser  in 
the  5-/X  region.  In  all  cases,  significant  scale-ups  will  be  required  to 
generate  enough  power  to  drive  a  full-scale  facility. 

Other  aspects  that  must  be  examined  include  optical  propagation 
through  the  high-temperature  air  and  the  possible  interaction  of 
high-power  laser  beams  with  the  walls.  Of  particular  concern  is 
clectric-field-generatcd  breakdown  in  the  energy-addition  region  of 
the  flow.  Although  this  occurs  at  very  high  intensities  for  laser  radi¬ 
ation,  the  variation  of  breakdown  as  a  function  of  pressure  is  of  con¬ 
cern.  It  is,  perhaps,  encouraging  to  note  that  the  breakdown  threshold 
reaches  a  minimum  with  pressure,  which  in  air  at  the  CO2  laser  wave¬ 
length  occurs  at  approximately  20-atm  pressure.  At  pressures  higher 
than  that,  the  breakdown  threshold  begins  to  increase  since  the 
breakdown  dynamics  become  dominated  by  collisional  processes.^ 
Breakdown  is  also  very  strongly  dependent  on  the  frequency  of  the 
electric  field,  with  the  breakdown  threshold  being  much  higher  for 
high  frequency,  i.c.,  near  infrared  and  visible  lasers,  than  for  low 
frequency,  i.c.,  far  infrared  lasers  or  microwave  devices.  Another 
complexity  that  will  have  to  be  investigated  is  the  interaction  of  the 
source  radiation  field  with  the  two-dimensional  flow  structure  in¬ 
cluding  the  boundary  layer.  The  high-tcmpcraturc  boundary  layer 
near  the  walls  will  cause  the  air  density  to  be  lower  there  than  in 
the  central  portion  of  the  flow.  This  density  variation  will  lead  to 
an  indcx-of-refraction  contour  that  may  be  used  with  the  geometry 
shown  in  Fig.  4  to  help  guide  the  lasers  and  avoid  excess  wall  heating. 

Alternate  methods  for  depositing  energy  downstream  of  the  throat 
need  to  be  examined.  A  particularly  attractive  approach  is  to  use 
very  high-power  microv^avc  gyrotrons  to  drive  molecular  oxygen  in 


the  vicinity  of  60  The  very  rapid  coupling  of  the  oxygen 

molecule  to  the  translational  kinetic  modes  in  air  suggests  that  this 
approach  can  potentially  be  used  to  deposit  large  amounts  of  energy. 
Microwave  gyrotron  sources  in  the  megawatt  region  arc  currently 
available,  and  an  increase  of  power  on  the  order  of  a  factor  of  1 0  can 
be  anticipated  in  the  near  future,  particularly  given  that  continuous 
operation  may  not  be  required  for  this  facility. 

Conclusions 

These  preliminary  examples  suggest  that  a  radiativcly  driven 
hypersonic  air  facility  may  be  feasible.  Major  questions  still  re¬ 
main,  particularly  regarding  the  effects  of  transport  phenomena  at 
high  pressure,  the  energy-addition  mechanism,  and  the  associated 
chemical  kinetics.  The  basic  concept  of  achieving  high  enthalpy  by 
using  ultrahigh  pressure  and  coupling  energy  in  downstream  of  the 
throat  does,  in  principle,  lead  to  a  wind-tunnel  apparatus  in  which 
the  air  remains  cool  compared  with  other  approaches.  High-pressure 
plena  have  been  constructed  in  Russia,  and  ultrahigh  pressures  have 
been  achieved  in  small-scale  laboratory  devices.  Issues  relating  to 
gas  transport  also  must  be  further  examined.  These  include  such  crit¬ 
ical  considerations  as  nozzle  erosion,  wall  cooling,  and  nonequilib¬ 
rium  phenomena.  Although  these  questions  will  require  significant 
research  to  be  done  in  a  variety  of  fields,  it  is  worthwhile  to  point  out 
that  the  development  of  a  true  air  hypersonic  test  facility  will  have 
a  major  impact  on  the  design  and  testing  of  hypersonic  vehicles. 
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ABSTRACT 

We  present  some  of  the  fundamental  spectroscopic 
and  dynamic  issues  relevant  to  radiative  coupling  of  large 
amounts  of  energy  into  supersonic  flow.  The  work  is 
motivated  by  the  Radiatively-Driven  Hypersonic  Wind 
Tunnel  (RDHWT)  concqrt,  which  has  been  suggested  as 
a  ^ssible  path  to  future  hypersonic  ground  test 
facilities.  The  RDHWT  relies  on  ultrahigh  plenum 
pressure,  followed  by  downstream  radiative  energy 
addition  to  achieve  high  Mach  number.  In  particular,  we 
focus  on  issues  of  molecular  absorption  and  collisional 
thermalization  in  high  pressure  air  at  moderate  (1.000- 
3,(XX)  K)  temperatures. 

Qumtitative  models  are  developed  to  predict  in&ared 
absorption  properties  of  high-pressure  air,  including 
spectral  line  mixing  and  collision-induced  absorption 
bands.  These  models  are  found  to  agree  well  with 
available^  experimental  data.  Microwave  absorption 
spectra,  including  oxygen  resonant  absorption  at  60 
GHz,  collision-induced  absorption,  and  absorpUon  by 
free  electrons  are  described,  and  a  simple  quantitative 
model  is  developed. 

The  key  pathways  for  vibrational  relaxation  of 
infrared  energy  are  identified,  potential  rate-controlling 
processes  are  revealed,  and  a  simple,  single  vibrational 
temperature  model  is  developed. 
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1.  INTRODUCTION 

The  ability  to  develop  hypersonic  technology  is 
constrained  by  the  lack  of  appropriate  ground  test 
facilities.  The  difficulty  is  easily  understood.  If  high- 
temperature  air  in  a  plenum  is  expanded  isentropically  to 
the  desired  Mach  number,  then,  to  produce  air  flow  at  a 
Mach  number  exceeding  12  at  a  static  pressure  and 
temperature  of  10’^  atm  and  3(X)  K,  respectively,  the 
plenum  pressure  must  exceed  4,0(X)  atm,  and  plenum 
temperature  must  exceed  6,(X)0  K.  At  this  high 
temperature,  air  cannot  be  contained,  and,  if  it  could,  it 
would  be  highly  dissociated.  For  this  reason,  run  times 
of  conventional  hypersonic  facilities  are  short  (of  the 
order  of  milliseconds),  and  the  flow  in  the  test  section  is 
highly  contaminated  with  NO  and  radicals. 

A  new  approach  aimed  at  circumventing  this 
fundamental  problem  has  recently  been  suggested  at 
Princeton.*  The  new  concept,  called  a  Radiatively- 
Driven  Hypersonic  Wind  Tunnel  (RDHWT),  relies  on 
ultrahigh  pressures  (20,(X)0-50,(X)0  atm)  and  moderate 
tonperatures  (900-1500  K)  in  the  plenum.  This  dense 
air  is  initially  expanded  through  a  nozzle  to  4,000-7,000 
atm  and  a  Mach  number,  of  say,  2.  Since  the  speed-of- 
wund  in  this  dense  fluid  is  much  higher  than  that  of  an 
ideal  gas  with  the  same  temperature,  the  kinetic  energy 
at  this  point  is  many  times  higher  than  it  would  be  for 
the  COTresponding  expansion  of  an  ideal  gas.  Additional 
energy  now  must  be  supplied  to  the  air  to  achieve  the 
necessary  total  enthalpy,  while  continuing  expansion  and 
acceleration.  The  heat  addition  should  also  proceed  at 
temperatures  such  that  the  entropy  reaches  that  of  the 
desired  test  section  conditions.  For  the  energy  addition, 
radiative  sources,  such  as  lasers  or  microwaves,  were 
suggested  in  Ref.  [1]. 

By  keeping  static  temperature  relatively  low  at  all 
points  of  the  flow,  the  RDHWT  should  be  capable  of 
significantly  decreasing  throat  erosion  and  the  production 
of  atypical  chemical  species,  and  allow  relatively  long 
run  times. 
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While  the  preliminary  modeling  results  presented  in 
Ref.  [1]  were  encouraging,  there  are  many  major  issues 
which  have  to  be  considered. 

In  the  present  paper,  we  focus  on  some  of  the  issues 
relevant  to  the  physics  of  energy  addition.  We  discuss 
fundamental  mechanisms  of  infrared  and  microwave 
absorption  in  high-pressure  air,  and  develop  quantitative 
models  which  can  be  compared  with  the  experiments. 
Pathways  of  molecular  energy  transfer  and  theimalization 
are  ctxisidered,  and  qualitative  effects  are  discussed. 

2*  IMTRARED  ABSORPTION  IN  HIGH- 
PRESSURE  ATR 

2.1.  Basic  mechanisms. 

Transmission  and  absorption  of  infrared  radiation  in 
the  atmosphere  has  been  extensively  studied  over  the  last 
few  decades.  Comprehensive  databases,  such  as  the  Air 
Forre  Geophysics  Laboratory  compilation,  are  readily 
available,  as  are  codes  for  line-by-line  high-resolution 
computation.2-3  However,  in  order  to  quantitatively 
evaluate  infrared  energy  addition  for  the  RDHWT, 
absorption  characteristics  need  to  be  known  for 
temperatures  up  to  a  few  thousand  Kelvin  and  pressures 
up  to  several  thousand  atmospheres.  The  existing 
standard  codes  cannot  handle  this  problem  adequately, 
because  they  do  not  account  for  physical  phenomena 
which  become  very  important  at  high  pressures.  Let  us 
t»iefly  consider  these  phenomena. 

Lids  mixing.  Widths  of  pressure-broadened 
vibrational-rotational  lines  (typically  0.05-0.08  cm- 
Vatm)  become  comparable  with  the  spacing  between 
adjacent  lines  (typically  2  cm  *)  at  pressures  of  the  order 
of  10  atm.  This  results  in  an  interference  between  lines, 
or  line  mixing,  of  rather  complex  nature.  To  better 
understand  it  qualitatively,  we  may  assume  that 
molecular  collisions  occur  at  random,  and  the  duration  of 
the  collision  is  infinitesimally  short.  Between  these 
brief  collisions,  a  molecule  is  emitting  or  absorbing  at  a 
certain  characteristic  frequency.  A  molecule  radiating  at 
a  given  frequency  before  a  collision  may  jump  to  a 
quantum  state,  which  radiates  at  a  different  frequency  as  a 
result  of  the  collision.  This  random  jumping,  or 
diffusion,  between  the  lines,  is  called  line  mixing.  The 
effect  strongly  depends  on  collisional  coupling  between 
rotational  quantum  levels  and  manifests  itself  in  the 
transfer  of  spectral  intensity  from  peripheral  regions  of 
the  vibrational  band  to  peak  regions  of  the  vibrational 
band.  What  is  observed  in  the  spectra  is  band  narrowing, 
with  higher  absorption  at  the  peaks  and  lower  at  the 
wings.  In  the  limit  of  very  high  pressure,  the  rotational 
band  structure  collapses  to  a  single  rotational  peak  for 
each  aUowed  spectral  branch  (P,Q,R,  etc.) 

effect,.  The  familiar  Lorentzian  lineshape  of 
pressure-broadened  spectral  lines  is  based  upon  the 
impact  approximation,  which  states  that  the  dephasing 
collision  occurs  instantaneously  fast.  This  assumption 


is  valid  when  Avx'C<l,  where  Av  is  the  deviation  from 
the  line  center,  and  x  is  the  collision  duration  (square 
root  of  optical  cross  section  divided  by  the  mean 
velocity).  Quantitatively  this  means  that,  even  for  an 
isolated  line,  the  Lorentzian  lineshape  is  valid  only 
within,  perhaps,  3-10  wavenumbers  from  the  line  center. 
Beyond  this  range,  the  lineshape  is  expected  to  be  sub^ 
Lorentzian.  This  effect  becomes  very  important  at  high 
pressures,  when  absorption  at  any  wavelength  within  a 
band  involves  significant  contribution  from  far  wings  of 
many  spectral  lines. 

Collision-induced  absorption  (TIAV  At  high 
pressures,  even  molecules  that  are  not  normally  infrared- 
active,  become  strong  absorbers.  To  understand  this 
effect,  known  since  the  1940s,  consider,  for  example, 
nitrogen.  Although  an  isolated  nitrogen  molecule  does 
not  possess  a  dipole  moment,  it  does  possess  a 
quadnipole  moment.  In  a  bimolecular  coUision,  a  dipole 
moment  is  induced  in  each  of  the  molecules  by  the 
quadnipole  field  of  the  other,  and  dipole  absorption 
becomes  possible.  Since  the  induced  dipole  moment  is 
modulated  by  molecular  vibration,  rotation,  and  relative 
translational  motion,  vibrational-rotational  (infrared), 
pure  rotational  (far-infrared),  and  translational 
(microwave)  absorption  bands  become  “allowed."  At 
moderate  piressure,  die  absorption  coefficient  increases  as 
density  squared.  A  typical  absorption  coefficient  at  the 
center  of  a  fundamental  band  is  of  order  lO'^cm'^Agt’^ 
At  very  high  densities,  however,  three-  and  four-particle 
collisions  become  significant,  partially  compensating 
the  two-particle  effect  and  somewhat  decreasing  the 
absorption  coefficient 

2.2.  Sub-Lorentzian  model. 

As  mentioned  above,  a  consistent  theoretical 
computation  of  absorption  spectra  with  line  mixing 
requires  knowledge  of  not  only  spectral  line  positions, 
strengths,  and  Iffoadening  coefficients,  but  also  of  matrix 
elements  of  the  relaxation  opierator  responsible  for 
rotational  coupling  between  quantum  states.^"® 
Formally,  exact  computations  involve  rather 
cumbersome,  albeit  standard,  inversion  of  large  matrices 
of  state-to-state  rate  data.  Since  the  rotational  relaxation 
rates  are  often  poorly  known,  this  formally  exact 
approach  is  still  subject  to  considerable  uncertainty.  An 
even  more  serious  problem  is  the  fact  that,  at  high 
pressures,  the  wing  effect  is  equally  as  important  as  the 
line  mixing.  Mathematically,  this  results  in  a  frequency 
dependent  relaxation  operator  which  presents  a  great 
difficulty  for  the  theory.'^-^  Despite  some  recent 
theoretical  developments,  the  theory  is  still  far  from 
being  able  to  provide  a  simple  and  accurate  predictive 
algorithm.'* 

At  the  same  time,  an  empirical  approach  has  been 
widely  used  for  inteipolation  and  extrapolation  of  infrared 
absorption  spectra.^' The  approach  is  based  on  a 
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semi-empirical,  sub-Lorentzian  lineshape,  with 
absorption  coefficient,  a,  at  frequency  v,  calculated  as: 


Sj  y, 

^  Y?+(v-v„.,f 


x(i 


(1) 


where  the  sum  is  taken  over  all  the  transitions  within 
the  band;  and  Sj ,  v^j  ,  and  are  the  line  intensity, 
posititHi,  and  Lormtzian  half-width,  respectively.  In  Eq. 
(1),  the  Lorentzian  lineshape  function  is  multiplied  by 
an  empirical  frequency-dependent  correction  factor,  x(lv- 
Vojl).  The  X  function,  which  accounts  for  line  mixing 
arid  wng  effects,  is  assumed  to  be  constant  for  all  lines 
within  a  given  band,  and  is  t)q)ically  approximated  by  an 
exponential  or  polynomial  expression  with  a  few 
adjustable  parameters.  The  parameters  are  evaluated  by 
least-squares  fitting  of  experimental  data 

To  represent  the  effect  of  spectral  intensity  transfer 
from  wing-to-peak,  the  x  function  should  smoothly  cut¬ 
off  with  increasing  Iv-VojI,  and  be  normalized  so  as  not 
to  ^ter  the  integrated  absorption  coefficient.  The  cut-off 
typically  starts  2-8  cm**  from  the  line  center,  dependent 
upon  the  Lorentzian  linewidth,  y^.  Recently,  a 
theoretical  justification  for  the  “corrected  Lorentzian”,  or 
X-function,  approach,  has  been  jx-ovided  by  Boissoles  et. 
al.,  and  by  Roney 

A  normalized  empirical  %  function  for  COj  bands, 
suggested  by  Fukabori  etal.,  is  given  by:^® 


evident  as  increased  absorption  at  line  center  and 
decreased  absorption  in  the  wings  as  compared  to  the 
Lorentzian  HITRAN  prediction.  Dotted  vertical  lines  on 
Fig.  1  correspond  to  some  of  the  output  frequencies  of 
the  HF  and  the  nuclear-reactor-pumped'  Xe  laser. 
Absorption  for  some  of  these  lines  will  clearly  be 
affected  by  band  narrowing.  The  sub-Lorentzian  code 
pa-edictions,  assuming  ti=  3  and  v„='2  cm'*,  are  indicated 
by  the  dotted  line  in  Fig.  1.  The  overall  agreement  is 
quite  good,  although  there  is  considerable  difference  at 
the  two  local  minima. 

In  order  to  extrapolate  to  higher  temperatures,  the 
temperature  dependence  of  q  and  v„  will  be  investigated 
in  future  experiments.  Theoretically,  line  mixing  and 
wing  effects  are  expected  to  be  less  important  as 
temperature  is  increased,  which  implies  an  inaease  of  v„ 
with  temperature.*^'^  Although  a  precise  theoretical 
analysis  has  not  been  made  from  the  physical 
interpretation  of  v„  as  dependent  upon  a  combination  of 
the  mean  rotational  frequency  and  the  reciprocal  of  a 
typical  duration  of  an  intermolecular  collision,  we  can 
expect  v„  to  be  roughly  proportional  to  the  mean 
rnolecu^ar  velocity,  and,  therefore,  to  scale  something 
like  T  .  For  the  present,  we  assume,  as  a  first-order 
approximation,  that  v„  scales  as  T'^I  Since  the 
absorption  is  less  inherently  sensitive  to  the  parameter 
q,  we  assume,  fw  convenience,  that  it  is  independent  of 
temperature.  This  will  be  investigated  experimentally  in 
the  near  future. 


X  =  C,  forlv-v„jl<v„ 


(v-vo.;r 


=  C  . 
K+yf 

forlv-Vo,il>v„ 


+yi 


I  Iv- Vo.il 


(2) 


where  Vn  and  q  are  empirical  parameters,  and  C  is  a 
normalization  constant  given  by: 


Fukabori,  et.al.  found  a  reasonable  agreement  with  their 
experimental  data  at  60  atm  of  pure  CO,  when  n=  3  and 
Vn,=  8  cm'*.l® 

We  have  adopted  this  empirical  lineshape  and  used  it 
in  a  simple  line-by-line,  high-resolution  transmission/ 
absorption  code.  The  new  code  will  be  calibrated  by 
least-squares  fits  to  experimental  absorption 
measurements  at  high  pressure  (up  to  1000  atm.)  and 
temperature  (up  to  1000  K).  These  measurements  are 
currently  in  progress.  A  preliminary  result  for  the  2.7 
micron  CO2  band  in  air  at  123  atm  and  room 
temi^rature  is  shown  in  Fig.  1,  along  with  HITRAN 
predictions.  From  Fig.  1,  the  band  narrowing  is  clearly 


2.3.  Collision-induced  absorption  model. 

The  theory  of  collision-induced  absorption  (CIA)  is 
quite  advanced.  An  extensive  review  is  given  in  a  recent 
monograph  by  Frommhold.^^  ^3j 

[14],  and  references  therein).  The  first  step  is  to 
calculate  the  vibration-rotation  line  center  positions  and 
strengths  nonempirically  using  quantum  mechanics. 
Translational  effects  are  accounted  foe  by  broadening  each 
transition  using  a  semi-empirical  lineshape  function, 
usually  with  four  adjustable  parameters.  The  parameters 
must  be  found  by  fitting  to  experiment.  For  low 
densities,  the  computations  are  more  or  less 
straightforward,  although  cumbersome.  The  biggest 
problem  in  the  theory  arises  at  high  density,  where  there 
is  significant  contribution  from  many-body  interactions. 
To  our  knowledge,  there  is  currently  no  established 
simple  algorithm  for  dealing  with  the  many-body 
problem  in  conjunction  with  line-by-line  spectral  theory. 

To  compute  collision-induced  absorption  in  the 
fundamental  band  of  nitrogen  over  a  wide  range  of 
pressures  and  temperatures,  we  have  developed  a  simple 
model,  which  is  essentially  the  theory  suggested  several 
years  ago  by  Joslin,  with  modifications  based  on  results 
obtained  earlier  by  Guillot  and  Birnbaum.^^’^^ 
Although  not  as  accurate  as  full  line-by-line 
computation,  the  model  is  physically  and 
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computationally  simple,  and,  as  we  will  see,  can  provide 
a  reasonably  good  agreement  with  experiments. 

Joslin’s  theory  uses  the  well-known  relation 
between  the  absorption  coefficient,  a,  and  the  spectral 
intensity, 

where  c  is  the  speed  of  light,  and  A  and  k  are  Planck’s 
and  Boltzmann  constants,  respectively.  Joslin’s  theory 
uses  classical  mechanics  together  with  the 
correspondence  principle  to  calculate  the  spectral 
intensity.  The  classical  expression  is  then  multiplied  by 
a  statistical  quantum  factor 

where  Av=v-Vfujj.  The  factor  S,  recoirunended  by 
Poll,  approximately  accounts  for  quantum  effects, 
enforces  microscopic  reversibility  (detailed  balance) 
requirements,  and  produces  significantly  asymmetric 
band  shape. 

The  spectral  intensity  is  the  Fourier  transform  of  the 
total  dipole  moment  (per  unit  volume)  autocoirelation 
function.  Since,  in  a  first  approximation,  vibrational, 
rotational,  and  translational  degrees-of-freedom  are 
uncoupled,  the  autocorrelation  function  can  be  written  as 
a  product  of  a  single-molecule,  vibration-rotation 
correlation  function,  and  a  two-molecule  translational 
term.  Consequently,  the  spectral  density  is  the 
convolution  of  vibration-rotation  and  translational 
^jectra; 


^(V)=  J^'IvibrDt(vX(v-V').  (6) 


The  vibration-rotation  spectrum  can  be  calculated 
exactly,  because  it  contains  only  single-molecule  terms. 
Tbe  analytical  expression  can  be  found  in  Ref.  [15].  To 
compute  the  two-molecule  translatitxial  spectrum,  Joslin 
suggested  to  apply  information  theory. The 
informadon  theory  provides  the  “least  biased”  spectrum, 
based  on  available  information  wi  the  spectral  moments. 
Thus,  the  dynamical  problem  of  calculating  auto¬ 
correlation  functions  is  reduced  to  the  calculation  of  the 
first  few  moments  of  /„(v),  which  is  done  by 
equilibrium  statistical  mechanics.  What  simplifies 
computations  further  is  the  fact  that  odd  spectral 
moments  vanish  in  the  classical  limit 

In  our  model,  we  made  two  modifications  to  the 
original  Joslin  theory.  First,  to  compute  the  spectral 
moments,  a  knowledge  of  the  two-particle  radial 
distribution  function  (rdf)  is  required.  When  density  is 


very  high,  on  the  order  of  liquid  density,  an  explicit 
fonnula  by  Goldman  can  be  used.^*  However,  that 
formula  is  not  valid  for  lower  densities.  We  used  the 
classical  expressicm  for  the  rdf  in  the  low-pressure  limit, 

rdf^exp[-i^),  (7) 

whoe  u(r)  is  the  Lennaid-Jones  pair  interaction  potential 
The  low-pressure  rdf  was  then  matched  to  Goldman’s 
formula  at  the  lower  limit  of  the  latter’s  validity. 

The  second  modification  regards  many-body 
encounters.  Even  for  the  zeroth  and  second  spectral 
moments,  an  accounting  should  be  made  for  three-body 
collisions,  which  reduces  the  accuracy,  since  the  three- 
particle  rdf  is  poorly  known.  For  higher  moments, 
encounters  of  four  and  more  particles  should  be  taken 
into  account,  which  makes  the  calculations  intractable. 
We,  therefore,  chose  to  ignore  fourth  and  higher 
moments  and  to  remove  Joslin’s  three-particle  terms  in 
the  zeroth.  Mg,  and  second,  M2,  moments.  Instead, 
we  substituted  a  simple  correction,  suggested  by  Guillot 
Md  Bimbaum.^^  The  two  moments  were  multiplied  by 
linearly  deaeasing  functions  of  density,  p.  Specifically, 
Mo  was  multiplied  by  (1-p/po),  and  M^-  by  (1- 
p/(2po)),  where  the  parameter  po  should  correspond  to  a 
density  somewhere  between  the  close-packing  and  the 
solid  density.  For  nitrogen,  is  8(X)-9(X)  Amagat. 

With  these  modifications,  we  used  Joslin’s  method 
with  the  molecular  parameters  (quadrupole  moment, 
polarizability,  moment  of  inertia,  etc.)  listed  in  Table  2 
of  Ref.  [15].  To  validate  the  model,  we  computed 
absorption  in  the  collision-induced  nitrogen  fundamental 
band.  The  nitrogen  fundamental  CIA  is  of  particular 
interest  because  this  band  overlaps  some  of  the  lines  of 
DF  laser.  The  results  are  shown  in  Fig.  2,  along  with 
experimental  data  of  Reddy  and  Cho.20  As  seen  from 
the  figure,  the  agreement  is  very  good,  especially 
considering  the  simple  nonempirical  nature  of  the  model. 

3-  MICROWAVE  ABSORPTION  IN  HIGH- 
PRESSURE  ATK 

3.1.  Basic  mechanisms. 

Assuming  air  above  50,000  feet  consists  of  79% 
nitrogen,  21%  oxygen,  and  has  small  concentrations  of 
CO2  (330  ppm)  and  H2O  (5-10  ppm),  there  are  three 
IHinciple  mechanisms  for  microwave  absorption. 

Oxygen  absorption  The  microwave  spectrum  of 
oxygen  is  due  to  a  series  of  magnetic  dipole  transitions 
within  the  E  electronic  ground  state.  Since  the  ^E  state 
has  an  unpaired  electronic  spin,  the  nuclear  rotation 
(rotational  quantum  number  JV)  couples  to  the  electronic 
spin  S=1  to  produce  a  total  angular  momentum  J  which 
may  take  the  values  /=N,  JV-1,  N+1.  Transitions  within 
the  triplets  form  a  broad  band  around  60  GHz 
(wavelength  X=5  mm).  Motivated  by  the  interest  in 
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millimeter  wave  propagation  through  the  atmosphere, 
this  band  has  been  extensively  studied  both  theoretically 
and  experimentally.  The  highest  pressure  at  which 
measuremoits  have  been  performed  is,  to  our  knowledge, 
51  atm  of  pure  oxygen  at  room  temperature  (Mingelgrin 
et.  al..  Ref.  [21]).  With  increasing  pressure,  the  band 
was  observed  to  broaden,  with  the  maxiiniim  shifting  to 
lower  frequencies.  As  a  result,  the  absorption  coefficient 
at  60  GHz  grows  linearly  with  pressure  up  to  about  10 
atm,  and  then  remains  nearly  constant  at  0.07-0.08  m'*. 
Theoretical  models,  developed  by  Gordcm,^  Mingelgrin 
and  Gordon, Rozenkranz,^'^  and  Smith,^  stress  the 
importance  of  line  interference  (mixing),  similar  to  that 
in  the  infrared  spectra. 

Collision-induced  absorntion.  This  effect  is 
^entially  the  same  as  that  discussed  in  section  2  for  the 
infrared  region.  Since  the  transient  dipole  induced  in  a 
molecule  by  the  quadrupnile  field  of  another  molecule  is 
modulated  by  the  relative  motion  of  the  molecules,  there 
is  an  absorption  directly  into  translation.  Both  nitrogen 
and  oxygen  exhibit  this  effect,  with  nitrogen  being  more 
important  due  to  its  higher  concentration  in  air.  There  is 
no  specific  resonant  frequency,  and  miaowave  radiation 
of  any  frequency  can  be  absorbed.  As  for  the  infrared 
QA,  the  absorption  strongly  increases  with  density,  and 
while  the  effect  is  negligible  for  normal  atmospheric 
conditions,  the  absorption  coefficient  in  liquid  nitrogen 
reaches  0.1-0.4  m‘*.26 

Eree-electron  absorption.  Even  if  a  relatively 
small  number  of  free  electrons  are  present  in  a  gas,  this 
may  lead  to  considerable  microwave  absorption.  The 
mechanism  of  an  electron  gaining  energy  from 
electromagnetic  wave  is  essentially  an  inverse 
bremsstrahlung.  This  also  may  be  thought  of  in  simple 
classical  terms  as  follows.  The  electric  field  of  the  wave 
imparts  a  directed  velocity  component  on  the  electrons, 
which  is  superimposed  on  their  random  thermal  motion. 
Collisional  scattering  of  electrons  by  molecules 
randomizes  this  drift  velocity,  adding  to  the  energy  of 
thermal  motion.  At  steady-state,  the  energy  flow  from 
the  field  to  the  electrons  is  balanced  by  electrons  losing 
energy  in  elastic  and  inelastic  collisions  with  molecules, 
resulting  in  gas  heating.  As  we  show  later,  fractional 
ionization  as  low  as  10’^  can  result  in  strong  absorption. 
This  weak  ionization  can  be  due  to  the  natural  thermal 
ionization  present  in  air  at  temperatures  of  about  2500 
K,  when  chemical  contamination  is  not  yet  very  high. 

Before  proceeding  with  quantitative  models,  it  might 
be  instructive  to  analyze  some  general  features  of 
microwave  absorption  as  compared  to  the  infrared. 
General  formula  for  the  absorption  coefficient  is  given 
by  Eq.  (4)  of  linear  response  theory: 


In  this  equation,  the  factor  (l-e”"'’^*'^),  proportional  to 
the  difference  in  population  of  lower  and  upper  quantum 
states  participating  in  radiative  transitions,  reflects  the 
net  balmce  between  abscxption  and  stimulated  emission. 
For  infrared  spectra,  at  low  or  moderate  temperatures, 

hv>  kT,  and  the  factor  (l-e"'""^'^^  is  of  the  order  of  1. 
In  contrast,  in  the  microwave  -region,  where  the 
frequency  is  3  orders  of  magnitude  lower  than  that  of  the 

infrared,  hv  «  kT,  and  (l-e"'*'"^''^)  =  hv/(kT)  «  1. 
Physically,  this  means  that  because  of  the  close  spacing 
between  states  participating  in  microwave  transitions, 
their  populations  are  almost  equal  even  at  room 
temperature,  and  absorption  is  nearly  balanced  by 
stimulated  emission.  The  general  formula  then  reduces 
to: 

hv  r,  , 

“•^v— /(v)  (8) 

For  molecular  absorption,  typical  values  of  microwave 
spectral  intensities,  /(v),  should  not  be  much  different 
from  that  in  the  infrared  bands.  At  high  enough 
pressures,  when  bands  are  very  broad,  /(v)  becomes 
approximately  constant.  Also,  at  constant  density, 
bands  broaden  with  temperature  approximately  as  T‘  , 
and,  therefore,  ^ectral  intensity  near  the  band  center 
decreases  as  .  Thus,  from  our  analysis  we  arrive  at 
three  general  conclusions: 

First,  because  microwave  frequencies  are  3  orders  of 
magnitude  lower  than  those  of  typical  infrared  bands, 
molecular  microwave  absorption  is  inherently  much 
weaker  than  infrared. 

Second,  at  high  pressure,  microwave  absorption 
coefficients  should  be  proportional  to  the  square  of  the 
frequeiKy. 

Third,  molecular  miaowave  absorption  coefficients 
should  scale  with  temperature  approximately  as  at 
constant  density,  resulting  in  significant  weakening  of 
absorption  with  increasing  temperature. 

These  conclusions  refer  only  to  molecular 
microwave  absorption.  Free-electron  absorption  behaves 
differently  because  of  an  inherently  different  behavior  of 
/(v).  For  example,  if  electrons  are  created  by  thermal 
ionization,  then  their  concentration,  and  consequently, 
their  absorption  coefficient,  exponentially  increases  with 
temperature.  In  addition,  efficiency  of  electron  energy 
gain  from  an  electromagnetic  field  depends  on  the  ratio 
of  the  field  frequency  to  the  electron-molecule  collision 
frequency.  When  the  period  of  electron  oscillation 
induced  by  an  alternating  field  is  much  smaUer  than  the 
time  between  collisions,  the  collisions  cannot  efficiently 
convert  the  velocity  of  the  fast  oscillations  into  random 
thermal  velocity.  For  this  reason,  efficiency  of  electron 
energy  gain  is  the  best  in  a  constant  electric  field  (zero- 
frequency  limit).  Therefore,  the  absorption  coefficient 
should  decrease  with  the  field  frequency,  and  free-electrcm 
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absorption  is  much  stronger  in  the  microwave  frequency 
range  than  in  the  infrared. 


3.2.  Quantitative  models  of  microwave 
absorption. 

Oxygen  absorption  As  mentioned  in  section  3.1, 
rigorous  theories  exist  to  predict  oxygen  microwave 
absorption.  However,  these  theories  are  very 
computationally  intensive.  In  addition,  they  were 
developed  for  moderate  pressures,  and  scaling  with 
temperature  was  not  analyzed.  At  high  pressure,  many- 
particle  interactions  may  be  important,  which  has  never 
been  considered.  Therefore,  in  this  paper,  we  use,  as 
first  approximation,  simple  extrapolation  of 
experimental  data.  As  mentioned  previously, 
experimental  data  of  Mingelgrin  et.  al.  has  demonstrated 
that  at  60  GHz  the  room  temperature  absorption 
coefficient  grows  almost  linearly  up  to  about  10  atm, 
and  then  saturates,  staying  nearly  constant.^!  Since  10 
atm  of  pure  oxygen  corresponds  approximately  to  50 
atm  of  air,  (and  working  with  density  instead  of 
pressure),  we  assume  a  linear  increase  of  absorption 
coefficient  with  density  up  to  50  Amagat,  and  constant 
absorption  coefficient  for  higher  densities.  The  value  for 
50  Amagat  of  air  at  room  temperature  was  taken  from 
experiment. 21  For  temperatures  other  than  room 
temperature,  scaling  (T/BOO)'*'^  was  applied. 

Collision-induced  absorption  Calculations  of 
coUisitm-induced  absoq>tion  coefficients  were  performed 
similarly  to  those  for  the  infrared  CIA,  desaibed  in 
section  2.3.  In  fact,  the  microwave  CIA  computations 
are  simpler  than  those  for  the  infrared.  Since  the 
absorption  is  purely  translational,  only  in  Eq.  [6] 
needs  to  be  calculated.  The  first  step  in  the  model  is  to 
compute  the  zeroth,  Mq,  and  the  second,  A/j.  moments 
of  the  spectrum,  which  is  done  by  a  simple  integration 
of  the  two-molecule  rdf,  gir): 

Mo  =  (l  -  )j^  (9) 

(10) 


where 


f  M, 


exp\ 


2J 


2jt^M, 


M, 


0_v2 


(12) 


“  =  -3^^o 


Mn 


2jtM. 


exp\ 


2  J 


2^ 

M. 


(13) 


Results  of  the  calculation  are  shown  in  Fig.  3,  where  the 
collision-induced  absorption  coefficient  is  plotted  versus 
density  at  various  temperatures  and  two  frequencies~60 
and  110  GHz. 

In  fact,  the  absorption  coefficient  for  the  CIA  need 
not  be  calculated  at  every  value  of  temperature,  density 
and  frequency.  We  have  found  that,  at  fixed  temperature, 
the  absorption  coefficient  is  approximated  very  well  by 
the  difference  of  two  terms,  the  first  being  proportional 
to  p  ,  and  second  to  p^.  Scaling  of  the  room- 
temperature  calculations  as  reproduces  results  for 
higher  temperatures  with  very  good  accuracy.  As  per  Eq. 
[13],  the  frequency  dependence  is  simply Therefore, 
the  absorption  coefficient  can  be  written  as: 


a  =  [ap^ (14) 

where  coefficients  a  and  b  are  found  from  calculations  at 
room  temperature  and  60  GHz. 

Free-electron  absorption.  General  formula  for  the 
absorption  of  an  electromagnetic  wave  in  an  ionized  gas 
is  given  by  Raizer  as:^^ 


v„ 

a  = - ^ 

me  CO  +vi 


(15) 


whCTe  nj  is  the  electron  number  density,  a)=27rv  is  the 
angular  frequency  of  the  electromagnetic  wave,  Vn,  is  the 
electron-molecule  momentum-transfer  collision 
frequency,  e  and  m  are  electrcHi  charge  and  mass,  and  c  is 
the  speed  of  light.  The  collision  frequency,  proportional 
to  the  gas  density,  depends  on  electron  energy.  For  air 
this  dependence,  according  to  Kroll  and  Watson,  isr^S 


v„fe)  =  7.0  10’^ 


£  +  0.1 

£  +  4.94*’ 


sec'‘. 


(16) 


(11) 

In  these  formulae,  p  is  the  density,  po  is  a  parameter 
discussed  in  section  2.3,  Co  and  a.  are  molecular 
quadrupole  moment  and  polarizability,  respectively,  k  is 
Boltzmann's  constant,  and  nz  is  the  mass  of  the  nitrogen 
molecule.  The  spectral  intensity  and  the  absorption 
coefficient,  a,  are  then  calculated  very  simply: 


where  £  is  the  electron  energy  in  eV,  and  p  is  the  gas 
density  in  Amagat.  For  the  average  electron  energy  of 
0.5  eV,  and  expressing  p  in  terms  of  the  number  density, 
n,  in  cm'%  we  obtain:  Vn,=  2.87x1 0  */i.  For  the 
microwave  spectral  range,  at  densities  above  1  Amagat, 
the  field  frequency  is  much  smaller  than  the  collision 
frequency,  a>«v„,  and  the  absorption  coefficient  reduces 
to: 
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a  =  3.7-10®-^?^-  m’*  (17) 

n 

The  frw-electron  absoiption  coefficient,  therefore,  is 
proportional  to  the  degree  of  ionization  of  the  gas.  Even 

small  ionization,  such  as  10'^,  would  result  in  a 
significant  absorption.  This  small  ionization  degree  can 
be  provided  by  the  equilibrium  thermal  ionization 
existing  in  air  m  temperatures  2,500-3,000  K.  At  these 
temperatures,  air  contains  2-4%  NO,  which  is  primarily 
responsible  for  ionization  because  of  the  lower 
ionization  potential  of  NO  (9.26  eV)  compared  with  that 
of  nitrogen  or  oxygen.  In  the  range  2,500-3,000  K,  at 
modoate  pressures,  the  time  needed  to  establish  chemical 
equilibrium  in  air  is  less  than  the  typical  flow  time 
through  the  energy-addition  region  of  RDHWT.'  The 
degree  of  ionu^tion,  therefore,  can  be  computed  using 
standard  equilibrium  thermodynamic  codes  such  as 
ASTRA.  For  example,  at  a  pressure  P=10  atm  and 
temperature  T=2,500  K,  n/n  =  1.23  x  10''°,  while  at 
T=3,000Kn/n  =  7.92xl0^ 

Tbe  degree  of  ionization  can  also  be  evaluated  using 
Saha  equation:^^ 


h  ) 


(18) 


where  /i+  and  are  number  densities  of  positive  ions 
and  ionizable  species,  g+and  g,.  are  their  respective 
statistical  weights,  and  I  is  the  ionization  energy. 
Assuming  n^—n^.,  and  that  NO  is  the  ionizable  specie, 
an  approximate  formula  for  the  firee-electron  absorption 
^  be  obtained.  However,  this  formula  does  not  tairp 
into  account  the  existence  of  negative  ions,  and  may 
somewhat  overestimate  the  electron  number  density. 

As  discussed  above,  molecular  absorption  in  the 
microwave  spectral  range  is  inherently  weaker  than  in 
the  infrared.  Consequently,  while  infrared  radiation  can 
be  absorbed  in  a  single  path  of  moderate  length  (order 
meters),  complete  single-path  absorption  of  the 
microwave  energy  would  require  path  lengths  in  excess 
of  100  m.  For  this  reason,  microwave  energy  addition 
will  require  some  form  of  resonant  cavity  design  where 
an  intense  electroma^etic  standing  wave  could  build-up, 
being,  in  a  way,  equivalent  to  a  multiple-path  system. 

Another  critical  feature  of  the  use  of  miaowave 
energy  is  that  the  threshold  power  for  dielectric 
breakdown  of  air  is  orders  of  magnitude  lower  than  that 
for  infrared  radiation.^^  At  normal  atmospheric  density, 
the  threshold  for  microwave  breakdown  is  of  the  order  of 
1  MW/cm2.  Fortunately,  the  breakdown  threshold 
increases  as  the  square  of  the  desnity.27  Consequently, 
the  microwave  power  cannot  be  directed  from 
downstream,  as  assumed  in  Ref.  [1],  since  breakdown 
would  occur  in  the  low  density  region.  The  breakdown 
would  create  a  plasma  where  microwave  energy  would  be 


rapidly  absorbed,  stimulating  undesirable  plasmid 
chemical  reactions  and  preventing  the  microwave  beam 
fiom  passing  to  the  high-pressure  region. 

One  possible  way  to  accommodate  both  the 
multiple  cavity  path  requirement  and  restrictions 
imposed  by  the  breakdown  onset  is  to  utilize  a 
segmented  cavity  approach.'*^  Slot  reactances  on  the 
tui^l  walls  would  act  like  mirrors,  effectively  creating  a 
series  of  resonant  cavities  in  the  energy  addition  region. 
Since  the  cavities  would  be  highly  uncoupled  from  one 
another,  and  each  cavity  could  be  excited  independently, 
this  approach  offers  the  capability  of  tailoring  the  energy 
addition  profile.  The  field  would  be  blocked,  or 
significanUy  reduced,  in  order  to  prevent  the  onset  of 
breakdown. 

4.  MECHANISMS  OF  MOLECULAR 
ENERGY  TRANSFER  ANn 
THERMAI.TZATTON 

If  an  infrared  laser  is  used  as  an  energy  source  for  the 

RDHWT,  the  energy  is  deposited  initially  into 
vibr^onal  degrees  of  freedom  of  air  molecules.  To  be 
efficiently  used  for  flow  acceleration,  this  energy  must 
be  quickly  thermalized.  However,  vibration-translation 
(VT)  relaxation  is  well  known  to  be  a  relatively  slow 
process,  requiring  many  collisions.  If  the  energy  does 
not  thermalize  rapidly  enough,  it  will  not  be  converted 
into  kinetic  energy,  and  a  significant  nonequilibrium 
will  exist  at  the  test  section.  In  addition,  vibrational 
disequilibrium  would  trigger  processes  of  quasi-resonant 
vibration-vibration  (VV)  energy  exchange,  which  may 
strongly  overpopulate  highly  excited  vibrational  states, 
and  stimulate  undesirable  chemistry.^^’^O 

Rates  of  VT  relaxation  and  energy  transfer  among 
lower  vibrational  levels  of  nitrogen,  oxygen,  and  CO2 

are  known  quite  well.30.3l  RecenUy,  models  have  been 
developed  for  relaxation  involving  highly  excited 
states.^2  Both  the  models  and  data  are  valid  for  low  and 
nioderate  densities,  where  relaxation  is  determined  by 
binary  collisions,  and  the  mean  free  path  is  much  larger 
than  molecular  diameters.  In  the  RDHWT  energy 
addition  region,  however,  initial  air  densities  are 
typically  several  hundred  to  several  thousand  Amagats, 
which  is  of  the  order  of  water  or  liquid  nitrogen 
densities.  In  such  a  fluid,  the  distance  between 
neighboring  molecules  is  comparable  to  their  diameter, 
and  many-particle  interacticxis  are  quite  possible. 

Over  the  last  decades,  many  experimental  and 
theoretical  studies  of  vibrational  relaxation  in  liquids 
have  been  perfonned.^3-4l  Qjjg  qj-  ujQgf  important 
findings  of  these  studies  is  the  good  accuracy  of  the 
isolated  binary  collision  (IBC)  model,  which  states  that 
the  probability  of  relaxation  per  collision  in  the  liquid 
phase  is  close  to  that  in  the  gas  phase.  The  explanation 
for  this  is  simple.  In  order  for  many-particle  interactions 
and  correlations  to  be  important,  the  interval  between 
consecutive  collisions  should  be  less  than,  or 
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comp^able  to,  the  collision  duration.  Since  most 
collisions,  at  thermal  energies,  are  inefficient  to 
vibration^ly  relax,  a  “successful”  relaxation  collision 
requires  kinetic  energy  an  order  of  magnitude  greater  than 
kT.  Since  such  energetic  collisions  are  quite  rare,  the 
time  interval  between  them  turns  out  to  be  considerably 
longer  than  the  collision  duration.  Hence,  inter- 
coUisional  correlation  is  not  signiflcant. 

To  apply  the  IBC  model,  one  should  properly 
calculate  the  average  collision  frequency  in  terms  of  two- 
molecule  rdf.^^  Because  of  the  “excluded  volume” 
occupied  by  molecules,  the  mean  free  path  can  be 
expected  to  be  shorter,  and  the  collision  ffequency-- 
higher  than  that  calculated  by  extrapolation  of  low 
density,  gas-phase  formulae.  Therefore,  while  we  can 
approximate  relaxation  rates  as  being  proportional  to 
density,  we  must  keep  in  mind  that  this  is  a  conservative 
estimate  and  that  the  actual  high  pressure  relaxation  time 
will  be  shorter  than  the  extrapolated  value. 

Let  us  now  consider  the  principal  mechanisms  of 
the  relaxation  process  in  high-density  air  absorbing  HF 
laser  radiation  in  the  2.7  micron  spectral  range.  The 
process  starts  when  a  CO2  molecule  absorbs  an  HF  laser 
photon  resulting  in  simultaneous  excitation  of  the 
bending  (V2)  and  asymmetric  (V3)  stretch  modes: 

CO2+hv,„„-4CO2r02i;  (19) 

To  supply  the  necessary  enthalpy  to  the  flow  (typically 
5-10  MJ/kg),  and  assuming  a  CO2  concentration  of  330 
ppm,  each  CO2  molecule  must  absorb  on  the  order  of 
10  photons.  This  absorption  must  occur  during  the 
residence  toe  in  the  energy  addition  region  of  the 
tunnel,  which  is  typically  0.1-1  msec.  The  resulting 
average  pumping  rate  is  10-10®  photons/sec.  If,  as  is 
likely,  the  energy  addition  profile  is  nonuniform,  the 
rate  may  be  even  higher. 

The  CO2  bending  mode  relaxes  rapidly,  with  rate 
constant  1.6  x  lO'^^cm^s  at  1000  K,  and  even  faster  at 
2000  K.31  ConsequenUy,  the  bending  mode  thermalizes 
in  2  ns  at  10  Amagat,  and  0.2  ns  at  100  Amagat,  which 
is  an  order  of  magnitude  faster  than  the  average  pumping 
rate.  However,  if  the  energy  addition  rate  is  high  in  a 
region  wh^  density  is  of  order  10  Amagat  or  lower,  a 
saturation  is  possible,  resulting  in  air  becoming  almost 
transparent  to  the  laser  beam. 

Direct  VT  relaxation  of  the  CO2  asymmetric  stretch 
mode  is  slow.  Instead,  energy  can  be  transferred  to  the 
bending  mode: 

C02(001)-i-M->C02(030)+M  (20) 

with  subsequent  thermalization.  However  the  rate 
constant  for  this  process  is  an  order  of  magnitude  smaller 
than  that  of  near-resonant  energy  transfer  to  nitrogen:^  ^ 


The  rate  constant  for  this  process  is  2  x  10'*®  cm^/s  at 
1000  K,  and  10’*^  cm^/s  at  2000  K.31  Therefore,  as 
with  the  relaxation  of  the  bending  mode,  the  CO2-N2 
transfer  is  faster  than  the  pumping  rate  at  densities  above 
100  Amagat;  however,  at  densities  of  the  order  of  10 
Amagats  or  lower,  the  transfer  rate  may  become  equal  to 
the  pumping  rate,  implying  the  onset  of  saturation. 

Once  the  vibrational  energy  is  transferred  to 
nitrogen,  direct  VT  relaxation  is  slow,  with  rate  constant 
of  only  1.1  X  10'^’  cmVj  at  1000  K.31  More  efficient 
is  the  energy  transfer  to  oxygen: 

N2(v  =  1)  +  02(v  =  0)->N2(v  =  0)  +  02(v  =  1)  (22) 

The  rate  constant  for  this  process  is  6  x  10‘‘®  cm^/s  at 
1000  K  and  2  X  10'‘‘'cmVj  at2000  K.31  The  toe 
scale  of  this  process  is,  therefore,  0.3  psec  at  1000 
Amagat,  1000  K,  increasing  to  1  psec  at  10  Amagat, 
2000  K.  Assuming  a  pumping  rate  of  10®  photons/sec 
per  CO2  molecule,  and  330  ppm  of  CO2  in  air,  a 
nitrogen  molecule  would  receive  quanta  from  CO2  only 
once  every  24  psec.  The  nitrogen-oxygen  transfer  can, 
therefore,  remove  energy  i^idly  enough.  However,  once 
again,  if  the  energy  addition  rate  is  high  in  a  region 
where  density  is  less  than  the  order  of  10  Amagats,  or 
the  temi^rature  is  below  2000  K,  the  transfer  rate  may 
beccHne  insufficient  and  vibrational  nonequilibrium  could 
occur. 

The  final  stage  of  thermalization  is  oxygen  VT 
relaxation: 

O2(v  =  1)  +  M->O2(v  =  0)  +  M  (23) 

with  rate  constant  1.6  x  10’*^  cm^/s  at  1000  K  and  10' 
cm  /s  at  2000  The  corresponding  time  scale  is 
23  nsec  at  1000  Amagat,  1000  K,  and  37  nsec  at  10 
Amagat,  2000  K.  Consequently,  this  relaxation  is  much 
f^ter  than  nitrogen-oxygen  transfer,  and  is  not  a 
liimting  step  in  the  overall  relaxation  process. 

The  following  qualitative  conclusions,  relevant  to 
IR  lasers  as  the  radiative  energy  source,  can  be  drawn 
from  our  analysis.  First,  in  the  high-pressure  region  of 
the  flow  in  the  RDHWT,  relaxation  is  fast  enough  so 
that  no  saturation  or  vibrational  disequilibrium  is 
expected.  Second,  in  the  moderate-density  (about  10 
Amagat)  region,  particularly  at  temperatures  below  2000 
K,  Iwth  vibrational  disequilibrium  and  saturation  are 
possible.  Third,  the  over^l  relaxation  process  has  two 
possible  bottlenecks”:  CO2-N2  transfer  is  responsible 
for  the  saturation  effect  which  would  significantly  reduce 
absoiption,  while  in  the  absence  of  saturation  nitrogen- 
oxygen  energy  transfer  is  the  rate-limiting  step. 


CO2  (001)  +  Nj  (V  =  0)  CO2  (000)  +  N2  ( V  =  1)  (21) 
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5.  RELAXATION  MODEL 

The  observations  and  estimates  made  in  the  previous 
section  allow  us  to  use,  as  a  first  approximation,  the 
following  simple  model.  A  single  vibrational 
temperature,  Ty,  is  ascribed  to  air  molecules,  and  both 
the  average  vibrational  energy  per  molecule,  Cy.  and  the 
relaxation  equation  are  written  in  the  harmonic  oscillator 
approximation.  The  vibrational  relaxation  time  is 
assumed  to  be  determined  by  the  nitrogen-oxygen 
transfer,  so  that: 


The  number  of  photons  absorbed  per  second  by  CO2 
molecules  is: 


''a - f - 

CcnAfiVi^„  [1  -  exp(- jj 

where  P  is  the  radiation  power,  a  is  the  absorption 
coefflcient,  A  is  the  cross  section  area,  and  Cc  is  CO2 
molar  fraction.  We  require  that  be  less  than  Vt, 
which  implies  the  following  conditicxi: 


(24) 


e  = 


(25) 


where  0=3353  K  is  the  N2  characteristic  vibrational 
temperature,  [O^]  is  oxygen  number  density,  Eq  is  the 
equilibrium  value  of  Ey  at  Ty=  T,  W  is  the  rate  of 
energy  absorption  into  the  asymmetric  stretch  mode  of 
CO2,  and  kf/^  is  the  temperature-dependent  rate  ctmstant 
of  the  nitrogen-oxygen  transfer,  equal  to: 


^N-o  =  7. 03  X 10“’°  jcm^  /s 

An  estimate  of  the  saturation  effect  can  be  made  as 
follows.  The  fraction  of  nitrogen  molecules  in  the 

ground  vibrational  level  is  jj,  and  the 

fraction  of  CO2  molecules  in  the  first  exc 

V3  mode  is  xexp^~ 


,  where 


fraction  of  CO2  molecules  in  the  first  excited  level  of  the 

0, 

,  1 1 X  exD\  - 

'  V, 

0c  =  3407  K  is  the  characteristic  temperature  of  the  V3 
mode.  Therefore,  per  CO2  molecule,  the  rate  of  transfer 
to  nitrogen  is: 


-  l^c-N  (T)Cn/i  •  [1  -  exp{-  %)j 

[l-exp(-®5^)]-ejp(-®5^)  , 


(26) 


wbere: 

kc.N=4.2A-lO-^^T^ 

exp[l0-^7’(8.84-10-^r-207)+-^j  cmVj 


is  the  rate  constant  of  the  CO2-N2  transfer,  Q  is  the 
molar  fraction  of  nitrogen  (approximately  0.79),  and  n  is 
the  total  number  density  of  air  molecules.^®-^* 


P<^ ACcCf,n^ hv fl - exp{j- jl 

r  3  (28) 

f-^^v(-yEy)]exp{-yry)kc-N  ■ 

If  this  flow  condition  is  not  satisfied,  then  the 
absorption  process  is  subject  to  saturation. 

Modeling  results  presented  in  Ref.  [1]  indicated  that 
very  significant  energy  addition  occurred  in  a  region 
where  the  density  rapidly  dropped  from  about  100  to  1 
Amagat.  This  is  a  quite  general  phenomena  reflective  of 
the  inherent  loss  in  stagnation  pressure  associated  with 
energy  addition.  This  provides  guidance  in  choosing 
optimum  parameters  for  energy  addition.  Figure  4 
shows  the  relaxation  time  plotted  against  temperature  as 
a  function  of  density.  The  shaded  area  represents  the 
range  of  reciprocals  of  typical  pumping  rates.  For  the 
energy  to  thermalize  successfully,  the  relaxation  time  at 
temperatures  and  densities  corresponding  to  the  high 
absorption  rates  should  fall  below  the  shaded  area.  If  the 
temperature  corresponding  to  the  regions  of  significant 
pumping  is  too  low,  this  condition  may  not  be  satisfied, 
and  vitrational  disequilibrium  could  occur. 

CONCLUSIONS 

Quantitative  models  for  prediction  of  high  pressure 
air  absorption  coefficients  in  both  the  infrared  and 
miCTowave  spectral  regions,  up  to  high  pressures,  have 
been  developed.  The  models  agree  well  with  available 
experiment^  data,  and  can  be  further  validated  and 
corrected  by  comparison  with  future  experiments. 
Analysis  of  molecular  energy  transfer  processes  in  air 
under  typical  RDHWT  conditions  relevant  to  infrared 
energy  addition  has  been  performed.  The  analysis 
suggests  that  N2-O2  and  CO2-N2  vibration-vibration 
energy  transfer  processes  play  key  roles  in  thennalization 
of  absorbed  energy  and,  under  certain  conditions,  may 
create  “bottlenecks”  in  the  relaxation  process. 
Microwave  energy  addition  in  the  RDHWT  is  influenced 
by  relatively  small  absorption  coefficients  and 
limitations  imposed  by  the  onset  of  electric  breakdown. 
This  suggests  the  use  of  a  segmented  cavity  approach  to 
achieve  the  required  energy  addition. 
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Fig.  1. 


Absorption  spectrum  of  air  at  123  atm,  300  K,  in  2.7-micron  CO2  band.  Vertical  dashed  lines  indicate 

strong  laser  lines. 
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Fig.  2. 


Absoipuvity  in  the  mtrogen  fun^ental  collision-induced  absoqjdon  band.  Solid  Unes-present  theory. 
Lhamonds,  triangles,  and  squares— exjjerimenti  data. 


Absorption  Coefficient  vs.  Density  at  60  GHz  Absorption  Coefficient  vs.  Density  at  110  GHz 


Fig.  3. 


Theoretical  predictions  for  microwave  collision-induced  absorption  in  air. 
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